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Fluctuation-dissipation relations, i.e., the relation between two-time correlation and linear re- 
sponse functions, were successfully used to search for signs of equilibration and to identify effective 
temperatures in the non-equilibrium behavior of a number of macroscopic classical and quantum 
systems in contact with thermal baths. Among the most relevant cases in which the effective temper- 
atures thus defined were shown to have a thermodynamic meaning one finds the stationary dynamics 
of driven super-cooled liquids and vortex glasses, and the relaxation of glasses. Whether and un- 
der which conditions an effective thermal behavior can be found in quantum isolated many-body 
systems after a global quench is a question of current interest. We propose to study the possible 
emergence of thermal behavior long after the quench by studying fluctuation-dissipation relations 
in which (possibly time- or frequency-dependent) parameters replace the equilibrium temperature. 
If thermalization within the Gibbs ensemble eventually occurs these parameters should be constant 
and equal for all pairs of observables in "partial" or "mutual" equilibrium. We analyze these rela- 
tions in the paradigmatic quantum system, i.e., the quantum Ising chain, in the stationary regime 
after a quench of the transverse field. The lack of thermalization to a Gibbs ensemble becomes 
apparent within this approach. 
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I. INTRODUCTION 

Motivated by recent experimental advances in the field of cold atoms, the theoretical study of the non-equilibrium 
dynamics of isolated interacting many-body quantum systems is currently receiving increasing attention [l|-[j] . Among 
the several questions that have been addressed, a central one concerns the way in which a macroscopically large isolated 
system evolving with unitary quantum dynamics from a generic initial state approaches equilibrium. In this work we 
discuss how such a problem can be effectively addressed from a novel perspective, inspired by analogous studies of 
non-equilibrium classical systems. A brief account of this study can be found in Ref. @. 

In classical systems thermalization is usually justified by advocating a chaotic dynamics that should ensure ergodicity 
in phase space and thermalization in terms of a microcanonical ensemble [f|. This implies (for systems with short- 
range interactions) that a large subpart of a much larger system thermalizes to a canonical ensemble. This condition 
is typically satisfied for generic Hamiltonians, even though the time required to reach such an ergodic regime might 
increase with the volume of the sample and thus be extremely large. Ergodicity implies that the time average of an 
observable for a single realization of the system coincides with its ensemble average. Within the statistical ensemble 
one derives exact relations that govern the dynamics of the system, such as the fluctuation-dissipation theorem. The 
way in which equilibration and ergodicity are understood for quantum systems evolving in a unitary (isolated) manner, 
together with the more appropriate way to rationalize the coarse grained description of their properties in terms of 
statistical ensembles, has been debated for the last 80 years 0-flH- 

Our aim is to revisit here the problem of thermalization in isolated many-body quantum systems with tools developed 
for the study of the non-equilibrium behavior of classical glassy systems [l3 - [l4| . Our study builds upon a large number 
of papers published in recent years. Basic questions as to whether a stationary state is reached and how this state 
could be characterized have been addressed in a number of simple models, including the one-dimensional systems 
reviewed in Refs. [H [l6j. The first picture which emerged was the following: non-integrable systems are expected to 
reach a thermal stationary state characterized by a Gibbs distribution with a single temperature. Intcgrable systems, 
instead, are not expected to thermalize. However, their asymptotic stationary state should nonetheless be described 
by the so-called generalized Gibbs ensemble (GGE) in which each conserved quantity is characterized by a generally 
different effective temperature fl7M23| . On the other hand, other works [2^ - [35j have shown (or at least argued) that 
this scenario could actually be significantly richer. 

Indeed, it was suggested in Refs. [30j that depending on the system's parameters and the specific quantity under 
study a conventional Gibbs ensemble might effectively capture some relevant features of the non-equilibrium dynamics 
even in integrable systems. In particular, a combined analytic and numeric study of the transverse field Ising chain (30l . 
l3~i| suggested that observables that arc non-local with respect to the excitation of the Hamiltonian, display the same 
relaxation scales as if they were in equilibrium at a finite temperature, at least for small quenches [2lT |3C|. Instead, 
local quantities such as the transverse magnetization do not show thermal behavior with the notable exception of 
quenches to the critical point. In fact, compared to non-critical quenches, the critical point shows some remarkable 
properties [3l|, H^, H3] that can be attributed to the gapless spectrum and the linearity of the dispersion relation at 
low momenta. This is clearly seen in the scaling limit of one-dimensional models with these properties, for which 
the results of Conformal Field Theory (CFT) [36|| suggest the emergence of a unique scale which plays the role of 
a temperature, at least for certain (non-generic) initial states (23l. l38j] . Accordingly, in what follows, we will focus 
specifically on critical quenches. Though it has been rece ntly shown that the dynamical properties of the transverse 
field Ising chain are eventually described by the GGE [2l|, [23, our aim here is to highlight this non-thermal behavior 
within a different approach, that can be extended to cases for which analytical solutions are not available. 

Most of the previous studies of thermalization in quantum systems compare the correlation lengths, the coher- 
ence times, or the expectation values of particular time- independent quantities with their values in equilibrium and 
extract in this wa y e ffective thermodynamic parameters such as an effective temperature of the out of equilibrium 
system 0, [T3, HHHE l29M3ll [H, [H[. However, these requirements may be too restrictive and/or insufficient to in- 
vestigate equilibration issues in systems with complex dynamics. The first two criteria are restricted to exponential 
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relaxation [39( , whereas the last one ignores the dynamics of the system. 
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In Gibbs equilibrium the (times-dependent) correlation function between any two observables is linked to the linear 
response of one of these observables to a linear perturbation applied to the other one in a model-independent way. 
Indeed, while the functional form of the correlation and linear response may depend on the pair of observables used, 
they can be affected by the spectral density of the bath, and they may of course be model-dependent, the relation 
between them remains unaltered and just determined by the temperature of the environment. This relation is the 
statement of the fluctuation-dissipation theorem (FDT) that involves only one parameter, i.e., the temperature of the 
system (for simplicity we assume that the number of particles is fixed). Quite naturally, a test of Boltzmann-Gibbs 
equilibration then consists in determining the correlation and linear response of a chosen pair of observables and to 
verify whether FDT holds for them. 

The analysis of fluctuation-dissipation relations (FDRs), i.e., the relation between correlation and linear responses, 
in classical dissipative macroscopic systems out of equilibrium has revealed a very rich and somehow unexpected 
structure. For instance, the spatio-temporal relaxation in classical glassy systems (or even non-disordered coarsening 
systems [40l - |4^ |) is very different from an equilibrium one with, e.g., breakdown of stationarity (aging effects) and 
other peculiar features. Still, the FDRs show that the dynamics can be interpreted as taking place in different 
temporal regimes each of them in equilibrium at a different value of an effective temperature with good thermal 
properties [12i4l5j . Similar results were found in quantum dissipative glassy models of mean-field type [43j . 

The main purpose of this contribution is to propose the use of FDRs as possible tests of (at least partial) equilibration 
in isolated quantum systems. In particular, this can be done by "measuring" independently the two-time symmetric 
correlation C AB and at the linear response R AB (defined in more detail further below) of two generic quantities A 
and B. On the basis of C AB and R AB in the frequency domain one extracts a frequency- and observable-dependent 
inverse "temperature" /3 A ^{ui) through the fluctuation-dissipation relation: 

h ImR AB (u) = tanh ( ^/M ^ C ab (uj). (1) 

The quantity (i A ^ (uj) provides important information on the possible equilibration of the system and on the various 
time/energy scales within which partial equilibration might occur. 

Concretely, we apply this idea to test equilibration in an intcgrable quantum system, the quantum Ising chain, 
quenched to its critical point, for which it has been argued that at least some observables could equilibrate in the 
usual sense of having their static and dynamic properties determined by a single global temperature. More precisely, 
we compute independently the correlation function and linear response of several pairs of observables. We note that 
despite the existence of many studies of the dynamics after a quench of the transverse field [2l| - [23l [30l l3lL [3rj 
none of them discussed the behavior of the linear response functions. For each FDR we extract a parameter (actually 
a time- or frequency-dependent function) that with a definite abuse of language we call " effective temperatures" . The 
analysis of these quantities, especially whether they are constant over certain time or frequency regimes and whether 
they coincide for different observables, will inform us about the (non-) thermal character of the dynamics. 

Let us emphasize that the idea of using FDRs to investigate thermalization properties in non-equilibrium system 
is completely general. Here, for illustration purposes, we apply it to a specific problem - the Ising model in a 
transverse field - in order to demonstrate that equilibration does not occur in this case, in spite of some evidence for 
the contrary, mentioned above. We expect our approach to provide an efficient test of thermalization also for non- 
intcgrable quantum systems, even though the characterization of their real-time dynamics is a very hard problem, 
often limited to small system sizes and short time intervals. 

The paper is organized as follows. In Sec. [ Ill we review several definitions of effective temperatures proposed in the 
context of quantum quenches [l7], EM HM l2ll. I30L l3~il l3r| and classical and quantum dissipative glassy dynamics (]~2| — 
EH l43l IFlj . Section IIIII summarizes those features of the static and dynamic behavior of a quantum Ising chain that 
are relevant to our study. Section HVl illustrates our results on the FDRs for several observables: the local and global 
transverse magnetization and the order parameter. Finally, in Sec. [V] we summarize our findings, discussing their 
implications and some ideas for future investigations. As already mentioned, a preliminary account of some of our 
results appeared in Ref. 

II. EFFECTIVE TEMPERATURES 

The evolution of a quantum system with Hamiltonian H (T) is ruled by the unitary dynamics 

iH±m)) = H(r)\m), (2) 

where T is some control parameter and \ip(0)) — \i/jo) an arbitrary initial state of the system. The initial condition is 
often chosen to be the ground state of the Hamiltonian H(Tq) corresponding to a different value of the parameter T. 
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In this case, one usually refers to the evolution in Eq. ([2]) as resulting from a quantum quench, i.e., from a sudden 
change r — > T of the parameter of the Hamiltonian. Alternatively, one might consider the case in which \i/j(0)) is not 
a pure state, e.g., it is a mixed state corresponding to the canonical distribution with Hamiltonian H (To) and inverse 
temperature (3o [52( • The initial state is then a generic excited state that does not correspond to an equilibrium state 
of the new Hamiltonian H (T) and right after the quench the system is in a non-stationary, non-equilibrium, regime. 

In order to investigate the possible emergence of an effective thermal behavior of the system after the quench, one 
can introduce effective temperatures on the basis of the behavior of various quantities. In Sec. Ill Al we discuss some of 
the possible definitions based on the (asymptotic) behavior of one-time quantities, such as those widely investigated 
so far in the literature on quantum quenches. In Sec. Ill Bt instead, we focus on the definitions based on two-time 
(dynamic) quantities that have been used in the context of glassy dynamics. We insist upon the fact that we still do 
not know whether the effective temperatures thus introduced can be attributed a thermodynamic meaning. 

A. Energy and constants of motion 

Among the various quantities that one can focus on in order to define an effective temperature, a special role is 
expected to be played by the energy of the system (II. l29j - [3ll . l33j : indeed the average energy E(t) = (ip(t)\H(T)\ip(t)) = 
(^o\H(T)\tpo) = E(t = 0) is conserved because the dynamics after the quench is unitary. Rather generally, one can 
define the density matrix p possibly describing the asymptotic state of the system as the one which maximizes the von 
Neumann entropy S[p] = — Tr[plog/5], subject to the constraint of having the correct expectation value of the energy. 
This amounts to assuming that the asymptotic state of the system long after the quench is effectively described by a 
Gibbs canonical distribution oc exp{— j3^ s H(T)}, in which the value of the effective temperature l//3f s is fixed by the 
constraint 

(MH(mo) = iTr[e-^ r )#(r)], (3) 

where Z is the partition function. The average on the l.h.s. is the energy of the system after the quench while the 
one on the r.h.s. is the average energy of an equilibrium state of H(T) at temperature T = T^ tl = l//3^ s . (Hereafter 
we set the Boltzmann constant ks = 1.) In this specific example, the l.h.s. of Eq. §3§ is independent of time because 
the energy is a constant of motion. 

In general, however, one would like to check that the temperature thus identified also describes the stationary limit 
of the average value of other observables. The time dependence of a generic observable O can be conveniently studied 
within the Heiscnberg representation 

6{t) = e idt de~ ikt , (4) 

within which the time-dependent expectation value on a generic mixed quantum state represented by a density matrix 
p (assumed to be normalized to one) is given by 

(6(t))=Tr[p6(t)}. (5) 

In canonical equilibrium at temperature T = /3 , p(T) is the Gibbs density matrix p(T) = cxp(— j3H) / Z(j3) and 
therefore, due to [H, p] = 0, the expectation value (0(t)) in Eq. is actually independent of time. In a generic 
non-equilibrium case, instead, the density matrix p over which the expectation value is taken describes the initial 
state of the system and it reduces to p = | , i/ , o)(' i / , o| when the system is initially prepared in a pure state \tpo)- In this 
case, analogously to what has been done above for the energy, one can compare the generic stationary values (if any) 
of the averages after the quench with the (time-independent) expectation value of the same observable O taken on an 
equilibrium canonical ensemble at the temperature T e ^ and then determine T e ^ in such a way that these two averages 
coincide, i.e., 

lim(TPo\6(t)\Tp ) = (6} T=T o (6) 

An effective thermal-like behavior of the system in the stationary state would require these temperatures TZ to be 
independent of O and, in particular, to coincide with Tjj defined above; however, this is not always the case |29l - [3l| . 

The discussion above implicitly assumes that the energy is the only quantity conserved by the dynamics and that 
in minimizing S[p] one has to account only for one constraint, which naturally leads to the Gibbs canonical ensemble. 
However, if the system is integrable, the situation turns out to be rather subtle because of the many "independent" 
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quantities which are conserved by the dynamics in addition to the energy [13; EE, IH| ■ Consider, for simplicity, a 
non-interacting Hamiltonian that can be written in the diagonal form 

H(T) = J2^(r)vlvk, (7) 

fc 

where the rjk's are creation operators for free (bosonic or fcrmionic) excitations of energy efc, labeled by a set k of 
quantum numbers. (For free theories k is the momentum.) The operators fik's satisfy the canonical commutation 
relations and the number hk of excitations is given by hk = f]]j)k- Clearly, [hk,H] = 0, i.e., the set {hk} is a set of 
constants of motion induced by -ff(r), independently of the initial state over which expectation values are calculated, 
which merely fixes the values of these constraints. Therefore, the dynamics after the quench are constrained by a 
large number of integrals of motion, the values of which have to be conserved. In repeating the minimization of S 
which leads to Eq. ([3]), it is necessary to introduce a number of Lagrange multipliers (see further below), one for 
each "conserved quantity", which eventually turn into a set of "effective temperatures" {T^} determined by the 
condition (V'ol^fclV'o) = ("fc) t=t\ { - These quantities prove to be particularly useful since they naturally appear in 
the calculation of (stationary and non-stationary) expectation values (T3, [H, ■ It was in fact suggested that 
the stationary behavior of the system after quenches towards Hamiltonians of the form can be described in terms 
of the density matrix oqqe obtained by maximizing the von Neumann entropy S[p] under the constraints on the 
expectation values of (hk)- This density matrix is of the form 

/3GGE = |e-^ A *'\ (8) 

where Afc = ek{T)/T^ s arc the Lagrange multipliers enforcing the values of the integrals of motion. 



B. Dynamic correlations and response functions 



As anticipated above, the aim of the present work is to introduce a definition of effective temperature which 
probes the dynamics of the system rather than the asymptotic time-independent properties discussed in Sec. Ill Al 
In particular, in this context, we are naturally led to consider FDRs, which turned out to be particularly useful 
in understanding various instances and features of the non-equilibrium dynamics of classical and quantum glassy 
systems [THl3|. 

The basic quantities which intervene in the FDRs are the two-time correlation between two generic operators A 
and B in the Heisenberg representation [see Eq. (UJ)] , defined by 

C AB (t, t') = (A(t)B(t')) = Tr[pA(t)B(t')]. (9) 

Clearly, for generic A and B one has (A(t)B(t')) ^ (B(t')A(t)) and it is natural to define symmetric and antisymmetric 
correlations as follows: 

C AB (t,t') = {[A(t),B(t')} ± ), (10) 

where [X, Y]± = (XY ± YX)/2. Without loss of generality we will consider either operators with zero average or we 
will imply that the average value is subtracted from the definition of the generic operator O: 0(t) —> (D(t) — (0(t)). 

In addition to C AB , another dynamic quantity of primarily importance is the instantaneous linear response function 
R AB which quantifies, up to the linear term, the variation of the expectation (A(t)) due to a perturbation which couples 
to the operator B, 



,ABU ^ - <W)> 



R (t, t') 



8h B {t') 



(11) 



where A(t) is obtained by evolving A - according to Eq. © - with the time-dependent perturbed Hamiltonian 
Hh B {t) = H — hs{t)B. In and out of equilibrium R AB (t,t') is related to the antisymmetric correlation C AB (t,t') 
defined in Eq. (fTOj) by the so-called Kubo formula [53j | . 

hR AB {t,t') = 2i6(t-t')C AB (t,t'), (12) 
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where 9(t) is the step function 9(t < 0) — and 9(t > 0) = 1 that enforces causality. In the following we will be 
concerned primarily with correlation functions of Hermitian operators, for which [C (t,t')]* = (B(t')A(t)) . Their 
symmetric and antisymmetric correlators can be expressed in terms of C AB in Eq. (|9|) as 

C AB (t,t')=ReC AB (t,t') and C A B {t 1 1') = i ImC AB (t, t 1 ), (13) 

so that Eq. QTZty yields 

hR AB (t,t') = -29(t-t')lmC AB {t,t'). (14) 

In equilibrium, the dynamics are invariant under time translations and therefore correlation and response functions 
are stationary, C± B (t,t') = C± B (t — t'), whereas out of equilibrium this is not necessarily the case. When dealing 
with stationary cases it is convenient to consider the Fourier transform of these quantities, for which we adopt the 
following convention: 

dt e iut f(t) and f(t) = / ^ e"*" */>)■ (15) 

-oo J — oo ^ 7r 

The canonical fluctuation-dissipation theorem (FDT) establishes a relation between the linear response of a system 
to an external perturbation and the spontaneous fluctuations occurring within the same system in thermal equilibrium 
at a temperature /3 . Remarkably, this relation does not depend on the particular system under consideration and 
takes the same functional form independently of the quantities which the correlation and the response refer to. 

In the case of canonical (Gibbs) equilibrium, the quantum "bosonic" FDT can be expressed in the time domain as 

R AB (t) = l -9(t) |^c—tanh C AB (u>), (16) 

where we reinstated h to make the classical limit h — > of the quantum FDT transparent. Indeed, in this case, one 
finds 

R AB (t) = -p6(t)±C AB (t), (17) 

which is the classical FDT. As expected, this limit is recovered for T = fi~ x 3> ft^typ where w typ is some typical 
energy scale of the quantum problem. The quantum FDT can be cast in a compact form in the frequency domain by 
Fourier transforming Eq. (jT6|) : 

~>AB i, \ _ . , ( fiAB { 



h Imi? AK (w) = tanh J C* b {uj). (18) 

Remarkably, knowing R AB and C AB for a pair of observables A and B allows the determination of the inverse 
temperature (3 of the system in equilibrium via Eqs. (|16[) and (|18[) . whatever the observables A and B are. In a sense, 
the FDT provides a viable method for "measuring" the temperature of a system, based on (local) measurements of 
correlations and response functions. 

Out of thermal equilibrium and in particular right after a quench in a isolated system the FDT is not expected to 
hold. It is tempting, however, to test whether FDRs such as Eqs. (ITBl and (|T8)) can be used to define a single "macro- 
scopic" temperature (or maybe a few), at least long after the quench and in the stationary regime. This approach 
turned out to be particul arly fruitful for understanding the physics of the thermalization of classical dissipative systems 
with slow dynamics [l3l 141. In particular, clarifying the relation between this temperature and the one defined from 
one-time observables 29h31| via Eq. ([3]) is definitely an important issue. In case some sort of thermalization occurs 
long after the quench, all these temperatures should become not only equal but also independent of the quantities 
used to define them. 

Depending on the specific quantum isolated system or model under consideration a stationary state may or may 
not be attained at long times. In several studies presented in the literature it was shown that a number of quantum 
isolated systems with short-range interactions reach a stationary state [l?], EH, H3, HH, Ull , while some fully-connected 
models [33L l35l . |54| and some mean-field approximations to models with short-range interactions [55l - [57| keep a non- 
stationary behavior. For the isolated quantum Ising chain a stationary state is indeed reached after the quench and 
we will therefore focus on this relatively simple case. Because of time-translational invariance of the stationary state 
we can equivalcntly consider two-time quantities in the time or in the frequency domain. According to the strategy 
outlined above, we define an effective inverse temperature (i AB (w) by enforcing the quantum FDR relation (|18[) . i.e., 

h ImR AB (oj) = tanh ( ^IfH h ^ (19) 
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where we consider C AB {t) and R AB (t) within the stationary regime. In complete generality f3 A g(u}) defined from 
Eq. (fill]) depends both on the particular choice of the quantities A and B which the correlation and the response 
function refer to and on the frequency w. Indeed, as the functional dependence of R AB {ui) and C ab {uj) on ui are, in 
principle, unrelated out of equilibrium it is necessary to allow for such a frequency dependence of f3 A ^ in Eq. (|19j) . 
The study of this dependence on lo provides an important piece of information on the dynamical scales of the system 
with respect to a given pair of observables A and B: heuristically, thermalization within a certain time scale would be 
indeed signaled by a f3 AB (uj) which becomes almost constant within the corresponding range of frequencies (T3, fl4|. 
This kind of analysis encompasses and generalizes in several respects previous studies in this direction. For example, 
in Ref. [58| an effective temperature was extracted for a system of one-dimensional bosons, after a quench of their 
interaction, by looking at the zero-frequency and zero-momentum limit of the FDR associated with the two-point 
density-density correlation function. Interestingly enough, such a temperature turned out to characterize the long- 
time and large-distance properties of the system after the adiabatic application of a periodic potential, for which 
a thermal-like behavior was found. In full generality, this analysis in the frequency domain can be extended and 
complemented by the analogous one in the time domain. Indeed, the time dependence of the response function R AB 
- which, in equilibrium, is connected to C ab (uj) via Eq. (|16j) — can be obtained from the Fourier transform of 
Eq. (fTT))) . Due to the integration over ui, the result of the possible variation of (3 AB with uj is "weighted" by the 
frequency dependence of C AB and therefore different " modes" contribute differently to the resulting time dependence 
of the response function. In order to highlight the possible emergence of "dominant" modes, one can define still 
another effective temperature B , based on the FDT in the time domain (p~6|) . i.e., 

R AB (t > 0) = * r dw e -fc,t tanh / Ki B fr^ 

hJ_ ao 7T V 2 / 

where, in contrast to Eq. (|19|) . the inverse effective temperature /3*g B on the r.h.s is assumed to be independent of 
uj. Note that Eqs. (fl"9| and (|20| are not equivalent, unless f3*^ B is allowed to depend on the frequency. In addition, 
while for given R AB (t) and C ab {uj) it is always possible to define the frequency-dependent inverse temperature 
(3 A ff(uj) from Eq. (fl"9|) . there might be no value of j3*g B for which the integral in the r.h.s. of Eq. (|20|) reproduces 
properly the functional form of the time dependence of the given response function R AB (t). Note that the frequency- 
dependent effective temperature f3 A g {uj) defined in Eq. (|19p is analogous to the mode-dependent effective temperature 
Pk introduced previously in the literature especially in connection with the generalized Gibbs ensemble. However, 
whereas the latter refers to the asymptotic averages of one time-quantities, the former accounts for the dynamical 
properties of the system in the stationary state. In the case of integrable models — as we will see below for the 
specific case of the isolated quantum Ising chain — these mode-dependent temperatures can be recovered as the 
frequency-dependent one /3^(u>) obtained by studying the correlation and response functions of specific observables. 

Within the stationary regime, it is rather natural to consider the behavior of the response and correlation functions 
at well-separated times. In fact, in classical coarsening and glassy systems this is the regime of structural relaxation 
in which partial equilibration of the slow (and non-equilibrium) degrees of freedom occurs fl3l Il4j . In consequence, 
we will focus on the effective temperature that emerges when one tries to relate the long-time stationary response 
and correlation functions after the quench via the long-time limit of the fluctuation-dissipation theorem in Eq. (|16[) . 
In particular, for large t the integral in Eq. (fH>|) is expected to be dominated by small values of w and therefore one 
can expand the hyperbolic tangent in a power series, which returns a sum over the odd time derivatives of C AB (t) 

nAB .. 2i ^ fiph\ n d n C AB (t) 1 d n tanh; 

R AB (t) = ~r J] c n 4- + where 



°° /AO fc\ " AnnAB I 



h <H V 2 / dt n nl dx r 

odd n— 1 



(21) 



x=0 



By inserting in this equation the expressions of the stationary response and correlation functions, R AB and C AB , 
after the quench one obtains an implicit definition of an effective temperature (3 = (3*g B which, hopefully, does 
not depend on time at the leading order and therefore provides a good definition of the temperature in the long- 
time limit. In Sec. IIV C 31 we will present an explicit determination of this temperature. At this point it is worth 
mentioning that if the correlation function on the r.h.s. of Eq. (|21[) decays as a power law at long times, then the 
leading order of the r.h.s. is indeed provided by the term with n = 1 and the possible temperature which one defines 
from this relation coincides with the one that one can infer from imposing (in the long-time limit) the validity of the 
fluctuation-dissipation theorem for a classical system, as in Eq. (|17p . However, as we will see in Sec. IIV[ oscillatory 
terms do actually modulate the algebraic decay of C AB ; accordingly, the leading term on the r.h.s. of Eq. (f2"T|) does not 
coincide with the first term of the expansion and the effective temperature (3 = (3*g B inferred from Eq. (|2 1 1) receives 
contributions from the derivatives of these oscillatory terms. Heuristically, the long-time behavior of the stationary 
response and correlation functions is expected to be determined by the low-frequency limit of their Fourier transform. 
In view of this fact, in the following we will be interested in understanding whether it is possible to recover the same 
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effective thermal description from the FDRs in the frequency and the time domains, at least in the low-frequency and 
long-time regimes. This is clearly possible in equilibrium where (3 is a constant. More precisely, in the non-equilibrium 
case we will compare the low- frequency limit of the effective temperature defined via Eq. (fl"9|) . i.e., 



2 , / HImR AB (uj)\ 

lim B^fP (u) = lim — arctanh — =- — 22) 

\ C ab {uj) J ' 



¥0+ e w->o+ Hlu 



with the value (3 = f3* A B obtained on the basis of Eq. (f2Tj) according to the procedure described thereafter. We 
mention here that in Sec. HVBI we will consider global quantities which are obtained by summing over all lattice sites 
and which are characterized by the fact that their correlation function C+ docs not vanish in the stationary regime 
even for well separated times. In these cases lm\ J ^ + C+(lu) ^ C+(oj = 0), and care has to be taken in evaluating 
the denominator of Eq. (|22l) . Finally, since one expects the quantum behavior to be relevant on the short-time scale 
whereas dccohcrcncc takes over at longer time differences, we will also consider the effective temperature extracted 
from the classical FDT in Eq. (fTT|). i.e., 



*B . Um _J_dC£(*) 

t™ BA B {t) At 



T AB = - lim nAR/ ^ ±^± , (23) 



with C AB and R AB taken in the stationary regime long after the quench. 

In passing, we mention that an effective temperature can also be defined on the basis of the relation between the 
linear response x (susceptibility) to a constant external perturbation Yia and the time- independent fluctuations in 
the thermodynamic conjugate quantity A [53 ], Indeed, in thermal equilibrium, the classical FDT in Eq. (|I7[) with 
A = B can be integrated in time, 



Sh A 



dt' R (t -t') = R aa (uj = 0) = j3C (t = 0) = f3[(A 2 ) - (A) 2 } (24) 



(where we assume C AA (t = — oo) = 0) which is indeed an alternative form of the classical equilibrium fluctuation- 
dissipation theorem. Note that, being the external perturbation Yla constant in time, the susceptibility x does not 
actually depend on the time t at which (A(t)) is measured. Out of equilibrium and in the stationary regime one can 
therefore introduce the additional effective temperature 

C AA {t = Q) , 

which has been used in Ref. [59| in order to test the thermalization of a specific isolated quantum system. Differently 
from the other quantities discussed so far, however, this effective temperature does not allow the investigation of the 
dynamic behavior of the system within different time- or frequency-regimes, as it involves only quantities which have 
been integrated out either in time or in frequency. 

III. THE ISING MODEL AND ITS DYNAMICS AFTER A QUENCH OF THE TRANSVERSE FIELD 

In this Section we briefly present the model, we recall its equilibrium phase diagram, and we discuss some of the 
properties of its dynamics after a quantum quench. 

A. The model 

We consider the quantum Ising chain in a transverse field jr > described by the Hamiltonian 

L 

where a°°' v ' z are the standard Pauli matrices acting on the i-th site of the chain, which commute at different sites. We 
assume periodic boundary conditions crf j+1 = af and we take the length L of the chain to be even. In what follows 
we set J, h, and kg — 1, i.e., we measure time in units of h/ J and temperature in units of JJks- It is well-known 
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that the Hamiltonian (f26|) can be diagonalizcd by performing three subsequent transformations. Firstly, we introduce 
Jordan- Wigner creation and annihilation fermionic operators ct . Cj [60l , l6lj which satisfy canonical anticommutation 
relations {cj, a,} = (5^ , {c$, cy} = {cj, ct} = and in terms of which 



The first identity implies = crt" + (cr^)^ and 



i=i 



(27) 



etc,. 

Having expressed all a x and tr z in terms of fermionic operators, the Hamiltonian becomes 

L 



H(T) 



L-l 

E 

i=l 



h.c 



44 



/i.e. 



(28) 



(29) 



where Np — 4^i 1S ^ ne number of fermions in the chain. The last term in this equation can be accounted 

for by extending the sum in the first term up to i — L after having defined cl+i = (— 1) Nf+1 ci, which amounts 
to assuming periodic boundary conditions for the chain if Np is odd and anti-periodic ones if Np is even. The 
Hamiltonian (|29[) conserves the parity of fermions and we restrict to the even sector which contains the ground state. 
Note that restricting to one of the two sectors is justified only when one considers expectation values of operators 
which arc defined in terms of products of an even number of fermionic operators, i.e., such that they do not change 
the parity of the state that they act on. 

The Hamiltonian H(T) in Eq. (|29[) , being quadratic, can be conveniently expressed after a Fourier transformation 



with k = ± 



?r(2n + 1) 
L 



and n = 0, . 



L 
' 2 



1. 



(30) 



where the sum runs over all allowed values of k. Note that in this expression we indicate both the fermionic operator 
c on the l.h.s. and its Fourier transform on the r.h.s. with the same symbol, the difference being made clear by the 
spatial (i, j) or momentum (k, I) indices and by the context. 

Finally, the Hamiltonian is diagonalizcd by a Bogoliubov rotation 



cos( 



Ck 
-k 



L -k 



where {7^} represent fermionic quasi-particles that satisfy the canonical anticommutation relations {7^,7^/} 
{7fei {'Ik') } = 3k,k', is a unitary rotation matrix and 



tan(20£) 



sin k 
r — cos k 



(31) 
and 

(32) 



For k > this relation has to be inverted with 26>£ <G [0, tt], whereas the values of 0\ for k < are obtained by using 
the property 9^ k = —0^. In terms of these quasi-particles the Hamiltonian H in Eq. (f29| reads 



£+(r) = 5>(r) (7^ + 7^-1 

fe>0 



where 



e fe (r) = 2Vr 2 - 2rcosfc + 



1 



(33) 



(34) 



is the dispersion law of the quasi-particles. In Eq. (|33[) the superscript + of H indicates that H + is the projection of 
the full Hamiltonian H in Eq. (|29|) onto the sector with an even number of fermions and that antiperiodic boundary 
conditions are enforced by choosing the wave- vectors k as in Eq. (|30l) . (Hereafter the superscript + is understood.) 
The ground state |0) r of the chain is the vacuum of the quasi-particles, defined by 7^ |0) r = 0, Vfc, which takes the 
form 



|0) r = 1] ( cos ^ + 1 sin ^ 4cl k ) |0) « J] e 

fc>0 



fe>0 



i(tan e k )cl&l k 



|6>, 



(35) 
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as a function of the fermions c k , where |5) is the vacuum of the fermions c k |5) = 0, Vfc. Hence, |0) r has the 

structure of a superposition of pairs c\,c*_ k , i.e., of pairs of fermions with opposite momenta. At zero temperature 
and in the thermodynamic limit, the system is characterized by a quantum phase transition at T = 1, where the 
gap of the dispersion relation e k (T) closes. The quantum phase transition separates a paramagnetic phase (PM, 
r > 1) with vanishing order parameter (<rf) from a ferromagnetic phase (FM, T < 1) with spontaneous symmetry 
breaking (of ) ^ and long-range order along the x direction. However, the long-range order disappears as soon 
as the temperature T takes non-vanishing values. As far as the transverse magnetization a z is concerned, instead, 
(of) ^ for all T > and all T > 0. 



B. Equilibrium and non-equilibrium dynamics 



Thanks to the transformations ([27]) . (|30|) and (|31|) the Hamiltonian H defined in Eq. (|26|) takes the quadratic 
diagonal form of Eq. (|33[) , which makes the model and its dynamics exactly solvable: indeed, in terms of the quasi- 
particle operators j k associated with H(T) one has access to all (thermo) dynamical properties. Within the Heisenberg 
picture these quasi-particles have a simple evolution 



= U{e k ,t)[ „' r fc t ), (36) 



[where e£ = £fc(r)] independently of the initial state of the system. Accordingly, the number operator h\ = 7^7^ 
of each kind of quasi-particle does not evolve in time n\(t) = n£ and its expectation value (n\(t)) on an arbitrary 
measure {e.g., on the initial state) is a constant of motion. 

In a quench, the system is prepared at t = in the ground state |0) ro of H(T ), while it is subsequently allowed to 

evolve, isolated, according to the Hamiltonian H(T). The quench from T to T injects into the system an extensive 
amount of energy which is henceforth conserved. (The statistics of this energy has been recently studied, e.g., in 
Refs. [ill l49j.) The dynamic observables one is interested in can typically be expressed in terms of the operators 
{&i}i,a and via Eqs. (|27[) and (|28|) in terms of the fermions {cfc,cL fc }fc in momentum space. A convenient way to 
calculate the associated dynamic correlations after the quench consists in expressing the (time-dependent) operators 
{ck(t),c i _ k (t)}k in terms of the operators {7^°}fc which diagonalize the original Hamiltonian H(To). The merit of 

this procedure is evident when calculating expectation values over |0) r , because {7^°}fc act trivially on their vacuum 
|0) Fo . On the other hand, the dynamics after the quench takes a particular simple form [see Eq. (|36|) ] if the operators 

one is interested in are expressed in terms of the quasi-particles of the final Hamiltonian H{T). Figure [1] summarizes 
schematically the relations between the various operators: black arrows indicate the transformations 1Z which connect 
them, given explicitly in Eq. (|31j) . The mapping in the direction opposite to the one indicated by an arrow is realized 
by the inverse transformation 1Zr l = . The grey vertical arrows indicate the time evolution, which takes the form 
of Eq. (|36|) in the basis of the quasi-particles {7^, 7^|.}. In order to solve the dynamics of the model, one first expresses 
these quasi-particles {7^, 7^} in terms of {7^°, 7_fc }j which requires a total rotation ^(6»[)^ t (^°) = TZ(S k (T,T )) of 
the suitable angle <5fc(r, r ) = 9 k — k °, as indicated in Fig.[T] Then the quasi-particles {7^, J_ k } are evolved according 
to Eq. (|3"6"j) in order to obtain {^(t), j_ k (t)}. In terms of the latter, the time-dependent operators {c k (t), c_ k (t)} arc 
eventually expressed according to Eq. (|3 X |) via a rotation 7Z^(9 k ). 

Combining these various transformations, the time-dependent operators {c k (t), ^L k (t)} are given in terms of 
{7^,7-^} by 



where 




Mel,t)K(5 k (T,T ))( * = * r Lfc \ J fc , (37) 

^(t) = e"^ 4 cos^ r cos(^ r - r k °) + e*E* sm6 r k sin(^ r - 6 r k °), 
<l' To (t) = ie- i£ fc t sin^cos(^ - 6£°) - ie 1 ^ 1 cos 6 T k sin(0£ - 6 r k °), 



and e£ = e k (T). This mapping allows one to express the average 



(•}=r o (0|.|0)r o 



(39) 
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FIG. 1: Schematic representation of the relations between the fermionic operators and the quasi-particles 7^ and 7^° of 
the final and initial Hamiltonians, respectively. Black arrows indicate the linear mapping provided by the Bogoliubov rotation 
Eq. (|31|l , whereas the grey vertical arrows indicate the time evolution, which takes the form (|36[) in the basis of the quasi-particles 
7^ of the final Hamiltonian -ff(F). 



over the initial condition |0) r in terms of ujv r °(i) and v^.' r °(t) defined above. 

After the quantum quench, all the obscrvablcs but the inte gra ls of motion (and possible functions of them) show 
a non-stationary behavior. After a transient (studied in Refs. j44l . |45| for the chain with free boundaries) the system 
reaches an asymptotic stationary regime. The typical time scale of this transient depends on the observable under 
study and on the initial and final values of the parameters, i.e., on Tq and T, respectively. For some observables the 
approach to the stationary value occurs via an algebraic decay in time. For other observables, instead, this decay 
is exponential and becomes faster upon increasing the energy injected in the system, which actually increases upon 
increasing |T — To|. 



C. Effective temperatures for the Ising model 



In this Section we specialize the various definitions of effective temperatures proposed in the literature for quantum 
quenches and discussed in Section III A[ in the case of the isolated quantum Ising chain. Because of the unitary 
dynamics, the energy of the system is conserved. Accordingly, if a thermal behavior emerges long after the quench, 
the (statistical) expectation value of the energy calculated on the corresponding ensemble has to match the (quantum- 
mechanical) expectation value of the energy of the system right after the quench [62|]. This suggests comparing the 
energy r o (0|-ff (r)|0)r o right after the quench with an equilibrium thermal average (H(T)) t=t e and defining an 

tE , 

LJJClctLUlC ± 

? .r 



effective temperature T o ^(r,r ) in such a way that these two averages coincide. By expressing H(T) in Eq. in 



terms of %° (see Fig. [TJ), one readily finds [2l|, |3fj, |3l| 

r 71 Ah 

r <o|#(r)|o>r =-J o ^ e fc (r)cosA,(r,r ), (40) 

where A fc (r,r ) = 24(r,r ) satisfies [see Eq. (J32J)] 

co,A t (r,r n ) = 4 ' rr °-' r + r °»7 t + " ,4!) 
e fe (r)e fc (r ) 

In Eq. (|40[) we took the thermodynamic limit L — > 00, which is assumed henceforth, and which allows one to 
replace J2k>o L Jo dk/(2ir). The angle Afc(r,To) in the previous equation is a crucial quantity, as it encodes the 
dependence on the initial state and fixes the (non-thermal) statistics of the excitations created at t = 0. Indeed, cos A& 
determines the expectation value (n£) over the initial state |0)r o °f the population n£ = 7^ 7^ of the quasi-particles 
of H(T) in the k-th mode 

, r, , p , 1 — COS At 

(nl) = (n r . k ) = k -, (42) 



which follows by direct calculation from Eq. (pTtj) . It is convenient to mention here that if the chain with Hamiltonian 
H(T) is in equilibrium within a Gibbs ensemble p(T) = exp[— H (T)/T]/Z(T) at temperature T, the average occupation 
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number (n£) = (fi^ k ) = l/[l + e efc ( r )/ T ] can be formally obtained from the expression (1421 (valid for the quench), 
with the substitution 

cos A k ^ tanh(e fe (r)/(2T)). (43) 

We will see below that this formal mapping is actually effective for a variety of observables, as it can be verified from 
direct calculations. In particular, the average energy within such an ensemble can be expressed as 

(H(T)) T = Tr [H(T)p(T)} = - f g efe (r)tanh (44) 

j o 

As in the case of the average occupation number, this expression can also be obtained from the corresponding one 
after a quench in Eq. (j40|) via the formal substitution in Eq. (|43]) . According to the strategy discussed in Sec. Ill Al 
one can implicitly define a global effective temperature Tj^(r,Fo) from the equality of energy averages 

r o (0|£(r)|0)r o = (£(r)) T=Teff(r>ro) , (45) 

where the l.h.s. is given by Eq. (|4T)|) and the r.h.s by Eq. (|4"4"|). 

In addition, it is also possible to define a mode-dependent effective temperature T^ s [13, HH Ho|, HH by requiring 
the integrands in Eqs. (|4"4"f and (|4T)]) to be equal, i.e., by defining T^ s such that 

cos A fc (r, T ) = tanh — . (46) 

v ' 2T c fc ff (r,r ) 1 ' 

This is nothing but the temperature that controls the population n\ of the fc-th mode and indeed T c fe ff could be 
cquivalently derived by imposing that each mode k is populated according to a Fermi distribution with temperature 
T^ S (T,T ) such that (n£) = (n^ fc ) = l/[l + e efc(r)/T *ff] . Note that in the thermodynamic limit T o fe ff (r,r ) becomes 

a continuous function of k. This means that the diagonal and quadratic structure of H(T) (i.e., the intcgrability of 
the model) naturally introduces an infinity (an extensive number) of "microscopic" temperatures, each one associated 
with a particular integral of motion The relevance of these temperatures is transparent by recalling that all 

expectation values are eventually determined by functions of However, for generic observables, these functions 

are typically combinations of multidimensional integrals, determinants, oscillatory factors in times, etc. and one 
cannot rule out a priori an emergent effective thermal behavior. 



D. Why a critical quench? 

Although all the considerations up to now are completely general, in what follows we will focus on the specific case 
of critical quenches, i.e., T = 1. This choice is motivated by the following heuristic arguments: 

1. There is some evidence that the expectation values of certain observables in the long-time limit after a quench 
to r — 1 do indeed coincide with the ones of a thermal state at temperature T^(T = l,To). This does not 
hold for quenches towards a phase with a gap (i.e., with r ^ 1). We will recall some of these results in the 
next Section. It is then natural to investigate up to which extent the apparent "thermalization" found for such 
one-time quantities at T = 1 carries over to two-time observables in the same stationary regime. The generalized 
FDRs precisely provide the tool to accomplish this goal. 

2. One can heuristically argue that the correlations between the entire system and a subpart of it might effectively 
act as a thermal bath for the subsystem. In this case, a gaplcss spectrum is expected to favor this effect, because 
energy exchanges at all scales are facilitated by the absence of a gap. 

3. Since our analysis will be primarily done as a function of the frequency u and eventually focus on the w — > + 
limit, it is natural to start the investigation in the absence of an energy gap in the spectrum. The gap might in 
fact determine a low-frequency threshold below which the spectral representation of some observables becomes 
trivial. Moreover, depending on the specific observable considered, the presence of the gap may introduce 
different frequency scales and non-analyticities in the Fourier transform of correlation and response functions. 
This issue definitely merits attention but it requires a dedicated study which is the natural continuation of the 
investigation presented here for T = 1. 
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FIG. 2: Quasi-particle occupation number (rife) as a function of k after the quench, for quenches towards the ferromagnetic phase 
F = 0.5 (left panel), the critical point V = 1 (central panel), and the paramagnetic phase V = 1.3 (right panel), and various 
initial conditions. We note that quenches within the same (gapped) phase are characterized by a low density of excitations hk 
at low energy, while quenches across the critical point are characterized by low energy modes with (rife) > 1/2 which correspond 
to negative effective temperatures T e *jj. 



4. The population (n£) of the k modes as a function of k [see Eqs. (|4"Tj) and (|42p] for quenches to (and from) the 
critical point is qualitatively different from the one involving the phases r, Tq =/= 1 with a gap in the spectrum, 
as shown in Fig. [5] Indeed, the distributions of (n fe ) for quenches at the critical point (or from it) have two 
properties in common with the equilibrium Fermi-Dirac distribution: 

• It is a monotonically decreasing function of the energy (i.e., of k for k > 0). 

• It varies within the range < (hk) < \ and therefore does not require the introduction of negative T e fe ff . 

When these conditions arc not satisfied (for quenches to the gapped phases), states exist at higher energy that 
have a larger overlap with the initial conditions than others at lower energy, and thus turn out to be more 
probable. This is in contrast with the behavior at the critical point and with any model at Gibbs equilibrium, 
where the probability of a given state is monotonically decreasing with its energy. 

5. The long-time, large-distance dynamical properties of a d-dimensional isolated quantum system after a quench 
from a ground state can be studied in terms of a suitable d + 1-dimensional problem in a slab (36^ . For a one- 
dimensional system quenched at a critical point with linear dispersion relation (i.e., dynamic exponent z = 1), 
this mapping has far-reaching consequences because the corresponding 1 + 1-dimensional problem is described 
by a boundary Conformal Field Theory (CFT) on the continuum. Interestingly enough, it turns out that the 
dynamic properties predicted within this approach are the same as those of the same CFT in equilibrium at 
a certain, finite temperature T, which should therefore naturally emerge after the quench of a generic one- 
dimensional system at its critical point. (Note, however, that the possible emergence of this thermal behavior 
depends on some rather general properties of the initial state [H[.) 



E. Dynamic observables 



In the following we will focus primarily on the transverse magnetization af and on the order parameter af , and we 
will denote by CJ(i) and C±(t) the corresponding autocorrelation functions. In addition, we will also consider the 

global magnetization M(t) = l/£Ei=i < j| an d the corresponding correlation function C±(t). These observables are 
distinguished by an important property [l], [3l| : af is non-local with respect to the quasi-particles in the sense that it 
has non-vanishing matrix elements with most of the states of the Hilbert space. On the contrary, af is local in the 
same variables, in the sense that it couples only few states. This distinction is more transparent if one recalls the 
expressions of these operators in terms of the Jordan- Wigncr fcrmions of Eq. (|27|) : of is a quadratic function of the 
Cfe and therefore also of the excitations 7^, while af is the product of a string of fermions and therefore it is non-local 
in the operators which diagonalize the Hamiltonian. 

Before presenting our results we briefly summarize what is known about the dynamics of C± x at the critical point 
T = 1 of the Ising chain in a transverse field. At equilibrium (r = T) the time decay of (af(t + tn)af (to)) as a function 
of t with fixed to is algebraic ~ |i| -3 / 2 at T = 0, whereas ~ at finite temperature [30]. In the isolated system 
after the quench (r ^ T), instead, the stationary decay of = §({of (t + t ) 7 af(t )}) — (af (t + t ))(af(t )) is still 
algebraic, but with a different exponent ~ \t\~ 2 [Hj]- When its initial condition is chosen to be fully polarized along 
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the ^-direction (corresponding to To = oo), C z _ — h([af(t + to),af(to)]) follows the same power-law decay ~ \t\~ 2 , as 
one can infer from the results of Ref. . 

The expectation value (&f(t)) of the order parameter of (t) decays to zero in the long-time limit for all r 7^ IV 
This is the same as in thermal equilibrium at T > 0, which is characterized by the absence of long-range order. The 
equal-time two-point correlation function (af (to)^ (to)} , instead, displays as a function of to an exponential relaxation 
towards its stationary value. This was first argued in Ref. }36j on the basis of semi-classical methods and Conformal 
Field Theory (CFT) and later shown to hold exactly via a suitable analysis of the model on the lattice [2l|. Moreover, 
\{af (t + to)af(to))\ was found [3(1 |3l| to decay exponentially cx exp{— t/rg} as a function of t (for fixed to), where 
the scale tq turns out to coincide numerically — at least for small quenches — with the one te which characterizes 
the time decay in equilibrium at a temperature T ~ T c g (with T c g defined according to Eq. (|35]l). The observed 
exponential relaxation is in contrast with the equilibrium correlations of the order parameter at T = 0, which decay 
algebraically both in space and time. To the best of our knowledge, instead, Cfi and therefore the response function 
R x have not been analyzed so far. 

IV. RESULTS 

In this Section we present and discuss our results about the behavior of a variety of correlation and linear response 
functions. This study allows us to assess the possible relevance of an effective thermal description of the dynamics 
following a critical quantum quench of the Ising model described in Sec. IIIII 

A. Transverse magnetization 

We start our analysis by considering the expectation value of & z (t) in the long-time stationary state after the quench 
(in view of translational invariance, we can drop the site index i from the notation; the case of a time-dependent T(t) 
was studied in full generality in Ref. [63^): 

dfc 

(a z ) Q = \im(a z (t)} = — cos(2^)cos A fc (r,r ), (47) 

t-K» J Q TT 

where here and in what follows the thermodynamic limit l/i^(i>o — * Jn dfc/ (2n) is taken. (This expression for (it 2 )q 
can be obtained straightforwardly by considering the time-independent contributions which emerge upon expressing 
the fermions in Eq. ([28]! in terms of 7^° via Eq. {37}, see also Refs. [H, [3J |H[.) As it is shown in Ref. [3l| this 
asymptotic value differs in general from the one (ct z }t that this observable would have in a Gibbs thermal ensemble 
at the temperature T = Tj| set by the energy of the initial state according to Eq. (|45j) . However, at the critical point, 
it turns out that [3l| 



' Q 



(48) 



which, by itself, would suggest an effective thermalization within a Gibbs ensemble at temperature Tj^. This is 
due to the fact that both for critical quenches in the stationary regime and at thermal equilibrium with r = 1, the 
expectation value of a z is related to the one of the Hamiltonian: (a z ) = — (H) /2. This relationship can be established 
by comparing Eq. (gO} and Eq. (gTJ, with the help of Eqs. (jAllj) and ([53]). which give e k (T) = 4cos(20£) for T = 1. 
Accordingly, Eq. (|48)) follows immediately from Eq. (|45|) . [In passing, we anticipate here that at the critical point, 
(a z ) Q = -E'(0), where E'(0) is given in Eq. (|A23]) .] 

The comparison between the expectation value obtained after the quench and the thermal one can be highlighted as 
follows (see also Fig. 9 of Ref. [3l[ for an equivalent analysis). For a fixed initial value Tq, Eq. (|4"5j) provides an implicit 
equation for T^(T, Tq) which can be solved as a function of T. The plot of this effective temperature is presented in 
the left panel of Fig. [3] as a function of T varying between and 3, for fixed r = 0.1, 1, and 2. As expected, the 
equilibrium value Tj^ = is recovered when T approaches the value of To corresponding to that particular curve. 
In addition, the curves clearly show that the effective temperature increases as the "distance" |T — To\ from the 
equilibrium condition increases. The knowledge of this effective temperature T ^(r,Fo) for a given pair of values 
(r,ro) allows one to calculate the expectation value that a generic observable would have at equilibrium in a Gibbs 
ensemble with temperature T = Tj^(r,ro). On the other hand, one can also determine the long-time expectation 
value (if any) of the same observable after a quench from To to T. In case thermalization occurs, these two values 
have to coincide upon varying T and To- On the right panel of Fig. [3] we present this test for a z , by plotting the 
thermal expectation value (ct z ) t=t e (vertical axis) versus the one after the quench (<j z )q (horizontal axis). The 
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FIG. 3: Left panel: Effective temperature T^(r,To) [see Eq. I|45p] as a function of T € [0,3], for fixed values of To = 0.1, 1, 
and 2. As expected, the equilibrium value Tjf — is recovered for r — ¥ To. Right panel: Comparison between the stationary 
expectation value (<5" z )q [see Eq. (|47[) ] after the quench and the thermal expectation value (<5" z )t in equilibrium at temperature 
T = T c f (r,r ). The solid lines are obtained by varying To 6 [0, 5], with fixed values of T = 0.5 (green), 0.75 (blue), 1 (red), 
1.25 (orange), 2 (purple), from bottom to top along the diagonal. The dashed lines, instead, correspond to fixed values of 
ro = (grey leftmost curve), 1 (red central), 5 (black rightmost) and are obtained by varying V. For critical quenches T = 1, 
(<t") t=t e and (<t z )q take the same value upon varying To and indeed the corresponding curve lies along the diagonal. 



various solid curves are obtained by varying Tq between and 5, with fixed values of T = 0.5 (green), 0.75 (blue), 1 
(red), 1.25 (orange), 2 (purple), from bottom to top along the diagonal. The dashed lines, instead, are obtained upon 
varying T between and 5, with fixed values of To = (grey leftmost curve), 1 (red central), 5 (black rightmost). For 
generic values of T and To the various curves clearly show that the average after the quench and the thermal average 
at temperature do not coincide, with the sole exception of quenches at the critical point r = 1, for which the 
corresponding solid curve in the right panel of Fig. [3] lies completely along the diagonal of the plot (thin black line) . 
As we anticipated in the previous Section, this result, apparently confirmed by an analogous observation in CFT [36| . 
was taken as an evidence for a sort of thermalization in the model, at least for this particular choice of parameters. 
However, it has been recently pointed out [23|, [38[ that the emergence of a well-defined, single effective temperature 
in CFT is actually a consequence of the particular choice of the initial condition in Ref. [36| . 

According to our general discussion in Sec. [TT] it is natural to investigate the issue of this apparent thermalization 
also in the light of the FDRs and of the various "effective temperature" parameters that can be extracted from them 
and which, in principle, depend on t or uj. These parameters are expected to provide a more sensitive tool to probe 
the possible asymptotic thermal behavior. 

The two-time symmetric connected correlation and linear response functions for generic T and Tq are given by (see 
App.[3]for the details of the calculation): 

/•*■ Ah f" At 

C* + (t + t ,t )=A — -Re[v k (t + t )v* k (t Q ) Ul (t + t )u*(t )}, (49) 
Jo t Jo n 

R z (t + t ,t ) = -8 6{t) — -lm[v k (t + t )v* k (t ) Ul (t + to)ut(t )], (50) 
Jo fJo " 

with Vk{t) = v^.' r °(t) and Uk(t) = u^' r "(t) given by Eq. (|38p . For critical quenches (F = 1, see App. IA 1[) and in the 
stationary regime to — > oo, one finds 

CUt) = lim CUt + t , to) = Jo(M) - E 2 (4t) + J*(U) - [E 1 ' {At)] 2 , (51) 

tQ— too 

R z (t) = lim R z (t + t ,t ) = 46(t)[J Ut)E(4t) - MAt)E' (At)], (52) 

to— »oo 

where here and in the following J a (r) indicates the Bessel function of the first kind and order a (see, e.g., chapter 10 
in Ref. [64|), whereas 

E(t) = / — sin(e fc r/4)cosA fc (l,r ), (53) 

J0 T 

[see Eq. (|A18l) and the plot of E(r) in Fig. Hp] with e k = e k {T = I) = 4sin(fc/2) [sec Eq. ([M|)]. The initial condition 
enters these expressions only via cosAfc [see Eq. (|41[) ]. Remarkably, at the critical point T = 1, this quantity turns 
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out to depend on To only through the ratio [see Eq. (|A12[) ] 

T=(}±£) 2 >1. (54) 

Note that T and consequently the correlation and response functions C\ and R z are invariant under the transformation 
To i — y T^ 1 which maps a paramagnetic initial condition into a ferromagnetic one and vice versa. However, this 
is true only for the stationary part of Ci and R z . Indeed, from the definition of 8k in Eq. (|32j) it follows that 

tan(26*^. + 26*£) = — tanfc and therefore, taking into account that Eq. (|32|) implies 26*£ =1 = (n — k)/2, one finds 

26l /T + 26\ = 7T - k, which yields A fc (r = 1, T ) = -A k (T = 1, T^ 1 ) and thus a change of the sign of sin A fc . With the 
help of the results in App. [3] especially Eq. (|A6[) . it is possible to see that the non-stationary terms in R z and Ci do 
actually depend on sin and, as a result, they are not invariant under the mapping Tq i— > . As we focus below 
only on the stationary regime, we can restrict our analysis to initial conditions in the ferromagnetic phase To < 1. 

In Fig.|3]we plot the stationary correlation [Eq. ([FT]) ] and linear response [Eq. ([52")) ] functions (left and right panel, 
respectively) as a function of time, in the case of the quench to the critical point T = 1, for To = (fully polarized 
case), 0.5 and 1 (equilibrium at T = 0). In both panels, the inset highlights in a double logarithmic scale the long-time 
algebraic decay of these functions (modulated by oscillatory terms), indicated by the thin dashed lines. If the system 
is initially prepared either deeply in the ferromagnetic phase T = or, equivalently, in the highly paramagnetic phase 
r = oo (both corresponding to T = 1), E(t) = Ji(r) (see App. I A 2| ) and Eqs. (jSTj) and (|52"j) become 



C%{t) = J 2 (4i) - -[J (At) - J 2 (4t)] 2 , (55) 



R z (t) = 2 0(t)J 1 (4t)[J o (4i) + J 2 (4i)]. (56) 

These results are consistent with the expressions reported in Ref. J3] for the symmetric correlation function C5_ and 
in Ref. [l6| for the response function R z after quenches starting from the fully polarized state To = oo, which are both 
generalized by our Eqs. (fSTj) and (|52|) . For a generic value of Tq (i.e., T ^ 1) E(t) cannot be expressed in terms of 
known special functions, though its asymptotic behavior in the long-time limit t 3> 1 can be determined analytically 
and yields [see App. IA 3l in particular Eqs. (|A37[) and (|A38[) ] : 



C z + (t) = -—co S (8t) + 0(r 3 ), (57) 



1 

8^rt 2 

R z (t) = t T_1 - sin ( 8t )] + °( t_3 )' ( 58 ) 

with T given in Eq. (|54|) . These expressions assume T _1 ^ 0, i.e., Tq ^ T = 1: indeed, in equilibrium To = T = 1 at 
zero temperature the relaxation is qualitatively different and the leading algebraic decay of both a nd R z [which, 
in this case, can be expressed in terms of Bessel and Struve special functions, sec Eqs. (|A30[) and (|A31|) ] turns out to 
be ~ t~ z / 2 , as discussed in App. IA 31 [sec, in particular, Eqs. (|A39|) and (|A40[) ] and highlighted in the inset of Fig. 2] 



In equilibrium at finite temperature, instead, such a decay is ~ t [31 1. 

In order to explore the possible definitions of effective temperatures based on the behavior of Ci and R z in the 
frequency domain, we consider the Fourier transform of Eqs. (|5T|) and (|52p . according to the definitions in Eq. (|15p . 
Due to the quadratic structure of these expressions in exp(±ie/ Cj ;) — in terms of which the trigonometric functions 
involved in the definitions of J a and E, see Eqs. (|A14p (|A19[) . are written — the corresponding Fourier transforms 
receive contributions only from real values of the frequency uj which coincide either with the sum + ei or with the 
difference — ti of the energies e^.i of two quasi-particlcs, depending on the range of uj. In turn, this structure is due 
to the fact that the observable a z under study is a quadratic form of the fermionic excitations {7^°, 7fc C ^}- Note that 
for |w| > w max = 2ek=n the Fourier transforms C+(lu) and R z (oj) of C^(t) and R z {t), respectively, do not receive any 
contributions from the integrals, because \ek ± ez| < 2efe =7r due to the existence of the upper bound at k = ir of the 
dispersion relation £fc(r). Accordingly, this results into a finite cut-off frequency ui max in the spectral representation 
of G z + and R z , and the corresponding Fourier transforms Ci_(uj) and R z (ui) vanish identically for |w| > w max . In the 
stationary state of the isolated system one expects the dynamics to be invariant under time reversal, which implies 
C±(i) = C±(—t) for the specific observable we are focussing on here; in turn, this implies the following symmetry 
properties 

C* b {uj) =C$ b (-cj) and Im R ab (lo) = -Im R AB (-u) (59) 

for the Fourier transforms of the symmetric correlation and response functions. In view of them, below we will restrict 
our analysis to the case uj > 0. In Fig. [S] (left panel) we present the result for Cl_(ui) (red) and R z (ui) (blue) for 
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FIG. 4: Time dependence of the correlation C+(t) (left panel) and linear response functions R z (t) (right panel) of the local 
transverse magnetization of [see Eqs. (|51[) and (|52[l ] in the stationary regime after a quench to the critical point T = 1. The 
blue, red and green solid lines correspond to To = 0, 0.5 and 1 (equilibrium at T = 0). The insets in the left and right 
panel highlight the algebraic decay of |C7+(£)| and \R z (t) — l/(47rT£ 2 )|, respectively, which are compatible with the asymptotic 
forms ()57[) and (|58l) for To 7^ T = 1 (blue and red solid lines), indicated by the lowermost thin dashed lines ~ t~ 2 in these 
double logarithmic plots. In the case of equilibrium at zero temperature ro = T = 1 (green solid line), both |C+(t)| and |i? z (£)| 
decay still algebraically, but with the different law ~ t~ 3 / 2 indicated by the upper thin dashed lines in the insets. See App. lA 3l 
and the main text for a detailed discussion of these functions and their asymptotic behavior. 



To = 0.3 (solid lines) and 0.8 (dashed lines), which we obtained by numerical integration of Eqs. (fSTj) and ([5^1> . These 
functions can be used to extract the frequency-dependent effective temperature T* s (oS), or equivalently /3* s (u)), from 
the FDR given in Eq. (flU)) . The dependence of the resulting T z s (uj) on the frequency ui is shown as solid lines in 
Fig. [5] (right panel) for various values of To, i.e., from top to bottom, To = (red), 0.3 (green), 0.5 (blue), and 0.8 
(yellow). The dashed horizontal lines indicate, in the same order from top to bottom, the corresponding (frequency- 
independent) values of the effective temperature Tjj = T^(T = 1,Tq) determined on the basis of Eq. (|45|) , as in 
Ref. [3l|, and in terms of which a thermal- like behavior was found for the asymptotic expectation value of of (see 
Fig. [3]). We point out that, beyond the fact that they seemingly tend to approach each other for uj ~ 4, as discussed 
further below, there is no obvious relationship between the effective temperature Tjj and the frequency-dependent 
T z s (uj), which even develops a mild concavity as a function of uj for 0.35 < r < 2.9. Note that T z s (uj) vanishes both 
for uj — > w max [with w max = 2efe =w (r = 1) = 8 at the critical point] and for uj — > 0, where it does so as 

Ufr, - 0) * 2 ' + (T-2)^an(v^T)/v^T (m) 

The comparison between this analytic expression and the actual behavior of T* s (uj — > 0) calculated numerically 
is shown in the inset of Fig. 1 of Ref. @. [Apart from this asymptotic behavior, C+(w), R z (uj) and T z s (uj) have 
been calculated numerically on the basis of Eqs. (|5ip and (|5"2"|) . Note, however, that these Fourier transform can be 
expressed as convolutions of those of J a (r) and E(t), the latter being discussed in App. IB 21 see Eqs. (|B29|) and 
(|B32|) .] The analysis of FDRs for the global transverse magnetization M(t), presented in the next Section, suggests 
that the vanishing of the temperature T z s (uj — > 0) can be traced back to the fact that the correlation and response 
functions for uj — > are not only determined by the low-fc modes, but they actually receive a contribution from the 
energy difference efc — ei between high-energy modes with k, I ~ n, which are indeed characterized by T fc f p 7r ~ (see 
the definition in Sec. IIII C|) . Accordingly, one can heuristically think of the behavior of C z {uj) and R z (oj) for u ~ 
and uj ~ 7r as being essentially determined by the presence of the lattice, which introduces a non-linear dispersion 
relation ek(T = 1) of the quasi-particles together with an upper bound 7r to the possible values of |fc|. In terms of this 
picture, the observation made above that T* s (oj) in Fig. [5] (right panel) approaches Tj| for u ~ 4 can be explained 
as a consequence of the fact that far from u> = and w = 8 more states contribute to the correlation and response 
function of a z which are less affected by the presence of the lattice, being in a sense closer to the continuum limit 
with linear dispersion relation, for which thermalization is expected (see point 5 in Sec. IIII D|) . 

For increasingly narrower quenches with Tq —¥ 1, T z s vanishes uniformly over all frequencies, as expected from the 
fact that the equilibrium value T — has to be recovered for T — >• V = 1. 
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FIG. 5: Left panel: Dependence of the correlation function (C+) and the imaginary part of the linear response function 
(Imi? z ) of <r z on the frequency uj, for To = 0.3 (solid lines) and 0.8 (dashed lines) at the critical point F = 1. C+(w) and 
Im R z (lj) are shown by a red (uppermost solid and dashed) and a blue (lowermost solid and dashed) line, respectively, and 
they vanish identically for |w| > w max = 8. [Note that the values of these functions for to < can be inferred from Eq. (1591) .] 
Right panel: Effective temperatures T* s (uj) = l//3 e ff(w) defined on the basis of Eq. (|19p . for T = 1 and various values of To. 
From top to bottom the red, green, blue and yellow solid curves correspond, respectively, to To = 0, 0.3, 0.5, and 0.8. The 
corresponding dashed horizontal lines indicate the values of the effective temperature T^g(T = 1, To) determined on the basis 
of the expectation value of the energy from Eq. ((45} . The comparison shows that there is no special relationship between 
these two possible effective temperatures, even though a thermal behavior was apparently observed when studying one-time 
quantities, i.e. with respect to Tj^ (see Fig. [3}- Note that for 0.35 < To < 2.9, T* g exhibits a mild concavity as a function of lu. 



We conclude that, although (df)g takes a thermal value [31|, the dynamics of of is not compatible with an 
equilibrium thermal behavior that would require all these temperatures to be equal within a Gibbs description. 
This case clearly demonstrates that assessing the emergence of a " thermal behavior" solely on the basis of one-time 
quantities might be misleading. 



B. Global transverse magnetization 

We consider here the global transverse magnetization 



(61) 



The generic two-time connected correlation and response functions of M(i) can be expressed as (see App. iBl) : 



v k (to)vk(t + t )u* k (t )u k (t. + t ) 



Cf (t + to,to)=8 C — Re 

JO 7T 

f n dfc 

R M (t + t 07 1 ) = -16 6{t) / — Im v* k (t )v k (t + t )u* k (t )u k (t + to) 

Jo 



(62) 
(63) 



where v k (t) = v^ r °(t) and u k (t) = u^' T °(t) are given by Eq. (|38|) . We point out that, in order to have a non- vanishing 
value when taking the thermodynamic limit of the corresponding expressions on the lattice, we have multiplied the 
fluctuations of M by a factor L, as it is required for global observables whose fluctuations are otherwise suppressed 
as L increases, see Eq. (|B3[) . For critical quenches (r = 1) the correlation and response functions in the stationary 
regime to — > oo read (see App. [B]) 



Cf(f)= lim Cf (t + to, to) = C+ MSt) + M8t) + F(8t) + F"(8t), 



t ->oo 



R M (t)= lim R M {t + t ,t ) = 46(t)[E(8t) +E"(8t)}, 



(64) 
(65) 
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where, in addition to the function E defined in Eq. (|53[) . we introduced the function F: 

r v Ah 

F(t)= — cos( efc T/4)cos 2 A fc (l,r ), (66) 

Jo T 

[see Eq_p7l) and the plot of F(t) in Fig. [TB] with e k = e£ =1 = 4sin(fc/2) and C = (1 + VT)- 2 [sec Eqs. (|BT0|) 
and (|B14[) ]. Interestingly enough, implementing the formal substitution of Eq. (H3J), i.e., cosA^. — ^ tanh(/3efc/2) 




FIG. 6: Time dependence of the correlation C+ (t) (left panel) and linear response functions R M (t) (right panel) of the local 
transverse magnetization M [see Eqs. (|64|) and (|65l) ] in the stationary regime after the quench to the critical point T = 1. The 
blue, red and green solid lines (top to bottom on the left panel, bottom to top on the right one) correspond to To = 0, 0.5 and 
1 (equilibrium at T = 0). The insets in the left and right panel highlight, in a double logarithmic scale, the algebraic decay of 
|C* 7 (t) — C\ and \R M (t)\, respectively, which are compatible with the asymptotic forms (|69p and (|70p for Fo 7^ T = 1 (blue and 
red solid lines) indicated by the thin dashed line in the left panel and by the lowermost thin dashed line in the right panel, both 
~ f~ 3/2 . In the case of equilibrium at zero temperature Fo = T = 1 (green solid lines), IC+^t)! still decays as t 3//2 whereas 
\R M (t)\ decays more slowly, as indicated by the uppermost thin dashed line ~ t 1 in the inset of the right panel. See App. IB II 
and the main text for a detailed discussion of these functions and their asymptotic behavior. 

in the definition (|B10[) of the constant C (see Eq. ([64)) ) docs not render its equilibrium value obtained with the 
equilibrium density matrix [rJ5| . This indicates that the dynamics of the model cannot be described solely in terms of 
the occupation numbers (see Sec. IIII CP — as the GGE does — because possible correlations (n£n£ fe ) between 
k and — k in the initial state can play an important role, at least for certain quantities P, [l8||. The constant C is 
in fact the result of time-independent correlations of the form c(ki,k2) = (^fci^fe 2 ) — (f l k 1 ){nk 2 ) Hi- While these 
terms do not vanish only if k\ = &2 within a statistical ensemble — such as Gibbs or GGE — which treats {fik}k 
as statistically independent variables, after the quench the summation on fci.2 of 0(^1,^2) which yield C receives 
contributions also from the term with ki = —k\, i.e., from c(fci, — k\) = c(k\, k\) due to the particular structure of the 
initial state. Generically, this latter contribution is subleading compared to the former as the system size increases 
and the thermodynamic limit is approached. However, this is not the case for M, due to its global nature. 

In Fig. [S] we plot the time dependence of the stationary (connected) correlation [Eq. (|6"4"]) ] and response [Eq. (p5|) ] 
functions of the global magnetization M (left and right panel, respectively) for a quench to the critical point T = 1, 
with Tq = (fully polarized case), 0.5 and 1 (equilibrium at zero temperature). In both panels, the inset highlights in 
a double logarithmic scale the long-time algebraic decay (modulated by oscillatory terms) of these functions, indicated 
by the thin dashed lines. 

As in the case of of discussed in Sec. lIV Al the initial condition enters the expressions of C¥ and R M in Eqs. (|64[) and 
([63)) only via the value of T defined in Eq. ([54"|) . In the relevant cases of an initial condition deep in the ferromagnetic 
or paramagnetic phase, T = or T = 00, one has T = 1, F(t) = [Jq(t) — Ja(r)]/2, and Eqs. (|M[) and ([63)) can be 
expressed completely in terms of Bessel functions [see Eqs. ([B18P and (|B19p and right before Eq. ([55|) ] 

Cf (*) = \ + |jb(8t) + lM8t) ij 4 (8t), (67) 
R M (t) = 0(f)[Ji(8i) + J 3 (8t)]. (68) 
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For generic To (7^ T = 1, i.e., T _1 7^ 0), the long-time decay of the stationary correlation and response of the 
global transverse magnetization for t > 1 are even slower than the ones of the corresponding quantities for the local 
transverse magnetization and is given by 

Cf(t) = C + sin(8* - tt/4) + 0(t~ 5 ^), (69) 

R M {t) = cos(8t - tt/4) + 0(r 5 / 2 ) (70) 

[see Eqs. (|B25p . (|B26[) . and (|B14|I ]. Note that the constant term C in Eqs. (|69|) and (|64| is not relevant for the purpose 
of studying the FDRs but it shows that the cluster property of the correlation function — which would require the 
connected correlation functions to vanish for well-separated times — does not hold because of the global nature of 
the quantity under study. The long-time behavior in Eqs. (|69[) and (|70[) is similar to the one observed at equilibrium 
at finite temperature [€351 ] . The case Fq = T = 1, i.e., T _1 = corresponds to the equilibrium dynamics at zero 
temperature T = for which the correlation and response function can be expressed as in Eqs. (|B21[) and (|B22[) in 
terms of known Struve and Bessel special functions, as discussed in App. [Bj In addition, as shown in App. IB 11 the 
leading algebraic decay of R changes into ~ f 1 , whereas the one of is the same as in Eq. (|69[) with C = 0. We 
also note that — differently from the case of of discussed in Sec. IIV Al — the leading-order decay ~ t~ 3 / 2 of both 
and R M for To 7^ T = 1 is actually the same as in equilibrium at T 7^ 0, which is discussed in Ref. [65J. Motivated 
by this observation, one could be tempted to extract the effective temperature of these dynamics by matching the 
features of the correlations and response functions after the quench with those of the same quantities in equilibrium 
at finite temperature, somehow extending to the present case the approach which was used in, e.g., the early studies 
of Refs. [3(1 l3l| . However, while the prefactor of the leading-order decay of both and R M in equilibrium depends 
upon T [65|, in the case of the quench the dependence on To appears only at the next-to-leading order, given that the 
long-i limit of the quantities in Eqs. and ([701) does not retain memory of the initial condition beyond the value 
of the constant C. 




CO CO 

FIG. 7: Left panel: Correlation function (C+ 1 ) and imaginary part of the linear response function (lmR M ) of the global 
magnetization M [see Eq. (|61[) ] as functions of the frequency ui, for To — 0.3 (solid lines) and 0.8 (dashed lines) at the critical 
point r = 1. C\' (ui) and lmR M (uj) are shown respectively by a red (uppermost solid and dashed) and a blue (lowermost solid 
and dashed) line. The dot-dashed black line shows the limiting value 1/2-^/1 — (uj/8) 2 of C+ and ImR AI for To — > 1. Due to 
the property stated in (|59|) we consider here only lo > 0. Right panel: Effective temperatures T^(uS) = 1/P^(oj) defined on 
the basis of Eq. (|19|) . for T = 1 and various initial conditions IV From top to bottom the red, green, blue and yellow solid 
curves correspond, respectively, to To = 0, 0.3, 0.5, and 0.8. The corresponding dashed horizontal lines indicate the values of 
the effective temperature T e g(r = 1, To) determined on the basis of the expectation value of the energy from Eq. (|45[) . As in 
the case of Fig. [5l the comparison shows that there is no special relationship between these two possible effective temperatures, 
even though a thermal behavior was apparently observed in the analysis of one-time quantities (see Fig. [3J). Note that for 
0.27 < To < 3.7, exhibits a mild non-monotonic behavior as a function of lu and that has a different qualitative 
behavior compared to T^ s reported in Fig. [S] 



In order to define a frequency-dependent effective temperature associated with M, we focus below on C¥ and R 
in the frequency domain. Differently from C+(w) and i? 2 (w), the Fourier transforms of Eqs. ([M)l and (p5[) receive a 
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contribution only from frequencies ui that equal ±2ek [in view of Eq. (|59|) we will restrict below to ui > 0]. This means 
that each frequency uj "selects" a mode k such that uj = 2ek (see App. IB 2l for additional details), and the high-energy 
modes with k ~ n do not contribute to the low-frequency behavior. These Fourier transforms reduce to integrals over 
k of terms of the form 5(uj — 2ek) which can be easily worked out: 

e fH^« + ivTHfi«?»»-H), (ri) 



Im R m (uj) = -^y/1 - (u/8) 2 - g(8-M)- ( 72 ) 

Note that, as in the case discussed in Sec. lIV Al these Fourier transforms vanish for \uj\ > uj max = 8. [For the purpose of 
defining the effective temperature according to Eq. (fT9]) we focus here — as we did in Sec. lIV Al — only on the imaginary 
part of R M (uj); its real part can be obtained from the Kramers-Kronig relation [67}, as discussed in App. IB 21 ] These 
expressions allow us to determine the effective temperature T^(uj) on the basis of the FDR in Eq. (fT9)) . Remarkably, 
this temperature turns out to coincide with the mode-dependent temperature TK — which characterizes the GGE 
in Eq. ([jo]) for the Ising model [U [23|, HH, [4(| — calculated for the value k u of k which is selected by uj, i.e., such 
that ui = 2efc^. Indeed, if in the definitions of E(t) and F(t) [see Eqs. ([55)1 and (|6"6"[) [ which appear in Eqs. ([6"4")l 
and (|65p one expresses cos A fc in terms of T e fe ff , the Fourier transform yields \LnxR M (w)]/Cf(ui) = tanh(w/(2T e fc ^)). 
Accordingly, comparing with Eq. (p~9|) . one concludes that T^(uj) = T e £ or, equivalently, 

T& = T%{u> = 2e h ). (73) 

Figure [7] (left panel) shows the frequency dependence of Ct(ui) (red) and Imi? M (u) given, respectively, by Eqs. (f7T|) 
and (J72J), for T = 0.3 (solid lines) and 0.8 (dashed lines). Note that Iim <4 ,_x)+ C¥ (uj) = 1/4 whereas lmR M (uj -> 0) = 
(VT"/16)w+C(w 2 ) and both functions vanish for a; — > ai ma x as (w max — w) 1 ^ 2 /4+C((w max — w) 3 / 2 ). This analysis in the 
frequency domain allows us to define the frequency-dependent effective temperature T^(uj) via the FDR in Eq. (|19l) . 
The function T^(uj) is shown as a solid line in Fig. [5] (right panel) for various values of To, i.e., from top to bottom, 
r = (red), 0.3 (green), 0.5 (blue), and 0.8 (yellow). The dashed horizontal lines indicate, in the same order from 
top to bottom, the corresponding (frequency-independent) values of the effective temperature = T^(T = l,Fo) 
determined as in Ref. [3l| on the basis of Eq. (|4"5j) . We point out that, as for T e 2 ff (w) reported in the right panel of 
Fig. [5j there is no obvious relationship between the effective temperature Tj| and the frequency-dependent T^(u>). 
Ttff(w) displays a non-monotonic behavior as a function of u for 0.27 < Tq < 3.7. Note that T^(uj) vanishes for 



UJ = UJ,, 



c as T* s (u>) does, whereas, at variance with it, T g(w) approaches a finite value 



-^eff = to T^)= * =2^1 (74) 

at low frequencies. [Note that taking the limit u> — > + is necessary in order to avoid the contribution cx S(uj) in 
C¥ (ui), see Eq. (|7ip .] This behavior is qualitatively different from the one of the frequency-dependent temperature 
T* s (oj) in Fig. [5J which vanishes for ui — > 0: the low-frequency value of T^g(u>) is solely determined by the low-energy 
modes which are characterized by a finite effective temperature T e fc ff . This property is unique of quenches to (and 
from) the critical point. 

Quite naturally, one may expect to recover this value by considering the FDR in the time domain for large 
times, as we discussed in Sec. lIIICl by replacing /? by a constant effective value (3* e in the r.h.s. of Eq. ([2"0"[) . the integral 
can be written as series of odd time derivatives of (t) , as in Eq. (|2 1 1) . Inserting the asymptotic long-time behavior 
of R M (t) and C¥(t) [see Eqs. ([6"5]) and ([70)1 ] in the r.h.s. and l.h.s. of this expression, respectively, yields the relation 
1 = tanh(4/3* ff ) at the leading order for t — >• oo and therefore T* s = l//3* ff = 0. The fact that T^g ^ T* s indicates that 
/3£fi(ui) = l/T^g(ui) cannot be approximated by an average constant in the integral on the r.h.s. of Eq. (|20l) . Indeed, 
since only the derivatives of the oscillating factor in Eq. ([6T?]) contribute to the leading order of Eq. (|20[) , T* e can be 
interpreted as a temperature associated with the oscillatory frequency uj = 8, corresponding to the threshold value 
w m ax- Consistently with this interpretation, the effective temperature T^(uj) defined on the basis of Eq. (|19p (see the 
right panel of Fig. [7]) vanishes upon approaching the threshold uj max , as T^(uj — > uj max ) ~ — 4/ln(w max — uj). The 
fact that this behavior turns out to be independent of T (and therefore of To) is also consistent with the fact that the 
leading-order term of the long-time asymptotic behavior of C¥ (t) and R M (t) [see Eqs. and ([70)) ] is — up to the 
constant C — independent of Yq and it is sensitive only to the largest frequencies (that, in turn, are associated with 
the largest energies). Such threshold w max results from the maximum of the dispersion relation and the quadratic 
dependence of M on the fermionic excitations, as noted for & z . From this analysis it is therefore unclear if one can 
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recover the value T^g = T^(lj — > 0) of Eq. (|74l) from an analysis of the response and correlation functions in the time 
domain. 

As we anticipated at the end of Sec. IIIB| an effective temperature can also be defined from the relation that 
connects in a classical equilibrium system the static susceptibility x to the fluctuations of the quantity which \ refers 
to, according to Eq. (|25|) (59| ; note, however, that this definition probes essentially the static behavior of the system, 
as the time dependence of the response function is "averaged" by the integration over time, which is known to miss 
the important separation of time scales in glassy systems. By using the results of App. iBland in particular Eqs. (JbTBJ) 
and (|B34|) . the effective temperature T c i )St for the global magnetization M can be expressed as 



rpM _C?(t = 0)-C?(t 

i cl,st _ 



-oo) 



R M (u) = 0) 



2VT(1 + VT) 



T - 1 



(75) 



Tife(l-T) - E e (l -T)' 



where K e and E e are the complete elliptic integrals of the first and second kind, respectively [see, e.g., chapter 19 in 
Ref. [H and the definitions in Eqs. (fA"2l| and (|A"25) ]. In Fig. [5] we compare this effective temperature (solid line) with 




FIG. 8: Comparison between the effective temperatures Tjj [see Eq. (|74[> . dot-dashed blue line] and T^f Bt [see Eq. (|75[) . red solid 
line] as functions of To or, equivalently, Tq 1 . For comparison, we report also the effective temperature Tjf (dashed black line) 
defined from the energy of the system. While all these temperatures vanish upon approaching the equilibrium case ro — ¥ T = 1, 
they are all generically different and in particular Tjf is always smaller than T~f st . 



T^fi [see Eq. (|7i)l . dot-dashed line] defined from the low- frequency limit of the frequency-dependent temperature T^(lu) 
(see Fig. [71 right panel), as functions of To (or, equivalently, r^ 1 ). While both of them vanish upon approaching the 
equilibrium case r — > T = 1 (i.e., T — > oo), they differ significantly for generic values of r , with T^f st approaching 
5/2 for To — > (equivalently, Tq — > oo, i.e., the fully polarized case T = 1), whereas — > 2 in the same limit. 
Accordingly, the temperature T^f st which relates the stationary fluctuations in the (stationary) state with the response 
to a constant external perturbation does not render the value T^g which is associated with the low-frequency limit of 
the dynamical properties of the system, and in particular it does not seem to have any clear relation with the various 
frequency-dependent temperatures T^g(oj) defined above for M. For the sake of comparison, in Fig. [8] we report also 
the temperature T^ s indicated by the dashed line. 



C. Order parameter 



In this Section we focus on the correlation function 

C x (t,t ) = {a?(t + t)af(t )) (76) 

[following the notation introduced in Eq. (f39|)] of the (local) order parameter of. The stationary correlation and 
response functions of the local and global transverse magnetization a z discussed in Sees. II V Al and ITV Bl turn out to 
be invariant — for quenches to the critical point T = 1 — under the mapping Tq \— > Tq . This is due to the fact that 
the dependence of the dynamics on Tq is brought about only by cos [sec, e.g., Eqs. (f5"Tj) . (|52l) . (|551) . (fMl) . (I6"5l) . and 
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(|66|) ] and cos Afc(l,ro) = cos Afc(l,rQ ) [see the discussion after Eq. (|54|) ] . In the stationary state attained long after 
the quench and for r = 1 we find numerically that this invariance also holds for C±. As a consequence, we will focus 
below on quenches originating from the ferromagnetic phase < Tq < 1. 



1. Computation of the correlation functions 

As it was already pointed out in Refs. [U [23|, [3(], |68|, the expectation value (a x (t)) of a x decays to zero at 
long times for any r ^ To. This is comparable to the equilibrium thermal behavior at finite temperature T > 0, 
which is characterized by the absence of long-range order of the magnetization along the x component. In order to 
verify the possible emergence of well-defined effective temperature (s) in the stationary state we focus on the two-time 
correlations CJ, which we computed on the basis of the method proposed in Refs. [13, [H| as an extension of a well- 
known approach in equilibrium [60j to the case of the dynamics after the quench. For periodic boundary conditions 
the computation of C± is non-trivial because the operator <rf (to + t)af(to) has non-zero matrix elements between 
states with different c-fermionic parity. Accordingly, for this observable, the assumption mentioned in Scc. lIII Al about 
the restriction to the even sector is not justified. However, following Refs. [3l|, Ho, HH], the correlation functions of 
this operator can be determined by computing a four-spin correlation function C xx , which can be done within the 
(antiperiodic) even sector, because C xx conserves the parity. On a chain of length L this correlation function is defined 
as follows: 



C xx (t, t ; L) = (o£ +1 (io + t) of (to) &l +1 (to + t) al {to)) 



(77) 



The spin cr£ +1 is identified with the spin a x after that the full string of Jordan- Wigner fermions, from 1 to L, has 
been inserted. By using the cluster property and by taking the thermodynamic limit, one can recover C±(t + to, to) 
from this quantity: 



[C x (t + t ,t )] = lim C xx (t,t ;L) 



(78) 



where C x (t 1 ,t 2 ) = C+(ti,i 2 ) +i C x (t 1 ,t 2 ). Following Ref. (SHH we introduce the operators A 3 (t) = c](t)+tj(t) and 
Bj(t) = ct(t) - Cj(t) in terms of the Jordan- Wigner fermions [see Eq. (|2"7])]. Note that (1 - 2c]cj) = AjBj = -BjAj 
and {Aj,Bi} = 0, Vj, /. Then, recalling the transformation in Eq. ([27]) we get 



C xx (t,t ; L) = ([&L +1 (to + t)---B L (t + t)][AL +2 (t + t)--- A L+1 (t 
x [B^to) ■ ■ ■ Bt(t )][A 2 {t ) ■ ■ • i i+1 (*o)]> 



(79) 



where the expectation value, here and in the following, is over the initial state, [Bi(t) ■ ■ ■ Bj(t)} = fj/ = i B lW) and the 
same holds for the product of operators Aj(t) . The next and final step is to write the Pfaffian in Eq. (f79")) in terms of 
the determinant of a 2L x 2L matrix [3lL kKJ, roll : 



[C xx {t,to;L)\ 



~{BA)\12 
-(BA)\l 



t +t,t 



( ba )£m 



(AA) 



-(AB)\l2 M (BB)%f a (BA)%£ 



ji 



-(AA)\ 0+t ' to 



■(BA)%% (AA) 1 ** 



(80) 



with (XY)* 1 / 2 = (X[i] (*i)Y[j] (t 2 )) and X,Y e {A, B}. The entries in Eq. §U\) are block matrices, whose indices 
j, I have subscripts which indicate their range, with ji,h & {■§ + 1, j + 2 . . . , L}, j 2 ,l 2 & {j + 2, \ + 3, . . . , L + 1}, 
J3) ^3 G {1: 2, . . . , -|}, and j'4, 14 € {2, 3, . . . , 4 + 1}. All the diagonal entries of the matrix in Eq. (|80| are zero, since 



24 



they do not enter the contractions. The matrix elements take the form: 

' {AjQtjMtt)) =^e^'- i )[ Ufc (t 1 ) + Ufe (f 1 )][<(t 2 ) + ^(i 2 )], 

k 

1 fc (81) 

(4(^)^2)) = I E eife °'~° ["*(*!) - U *(*l)] [^(*2) + ^(t 2 )], 

(4-(<i)A(*a)> = S^"" - w *(*i)] Kfe) - uj(*a)] , 

and the sums over k are consistent with the antiperiodic boundary conditions, i.e., they run over the values of k 
given in Eq. (|30[) . Some care is required in extracting numerically the real and the imaginary part of C x (t) from the 

value of [C xx (t,t ;L)] . In our case, starting form the initial value C x (0) = 1, we followed numerically the analytical 
solution in time. In passing, we mention that a different approach based on the solution of differential equations has 
been proposed in Rcf. [69| for computing similar correlation functions. 

2. Numerical results for the dynamics 

We computed C±(t + to, to) for an isolated chain with anti-periodic boundary conditions, L = 10 3 and to = 10. It 
turns out that these values of to and L are sufficiently large for accessing both the stationary and the thermodynamic 
limit (at least for values of To not too close to 1, see the discussion after Eq. (|87f further below), as we verified by 
comparing with results obtained for different values of <o an d L. Accordingly, we focus below only on the dependence 
on t o{C^(t + t ,t ). 

We analyzed the effect of the boundary conditions by studying also the dynamics of an open chain with free 
boundary conditions, a problem that has been thoroughly investigated in Refs. [HI, H^]. This case differs from the 
one considered here because in the open chain there is no distinction between the even and the odd sectors and 
translational invariance is absent. In spite of these differences, we found that, for fixed parameters, C± calculated in 
the periodic chain and in the bulk of the open chain perfectly coincide, at least for times such that finite size effects 
are not relevant. At larger times, instead, the correlation function manifests finite-size effects which are expected to 
depend on the boundary conditions; however, the study of these effects is beyond the purposes of this work. 




FIG. 9: Time dependence of the correlation C+(t) (left panel) and linear response R x {t) (right panel) functions of the order 
parameter, in the stationary state after a quench to the critical point Y = 1 and for different values of IV From bottom to 
top, red, green, blue and yellow lines correspond to Yq = 0, 0.3, 0.5, and 0.8, respectively. The insets highlight the short-time 
behavior of these functions and the dashed curves correspond to the analytic expressions in Eqs. (|82l) and (|B3p with Y — 1. 

In Fig. |wc show the time evolution of C$(t) (left panel) and R x {t) = 2i6(t)C x (t) [right panel, see Eq. (fHI) ] 
for various values of the initial transverse field Tq. (As anticipated, here and in the following we focus only on 
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quenches originating from the ferromagnetic phase Tq < 1 because the numerical data provide strong evidence that 



the stationary dynamics with T = 1 is invariant under the mapping r 



The insets in both panels provide 



a closer view of the initial-time regime and the corresponding dashed lines represent (with T 
behavior 

c x + {t) = i-2(n) 2 + o(t 4 ), 

R x (t) = A(a z ) Q Tt + 0{t 3 ), 



1) the short-time 

(82) 
(83) 

determined analytically in App.JCj with (cr z )g given in Eq. (|47|) . 

Qualitatively, the curves in Fig.[9]are intriguingly similar to those found for dissipative disorder spin models j43ll70l|. 
in the sense that they consist of a superposition of a monotonicl y d eca y ing tail and oscillations with decreasing 
amplitude. However, in the disorder spin model treated in Refs. [43l . |7C| there are the extra ingredients of aging 
phenomena and slow relaxation in the ordered phase, which are absent in the present problem. 

For t > 15, our numerical results are fitted very accurately by 



cut) 

R x {t) 



e- t/T A c 
e~ t/T A R 



ac 

t l/2 

an 

£1/2 



sin(4r, + </>) 
cos(4< + <f>) 



(84) 
(85) 



where Ac,r, ac,fl, t, and <fi are the parameters of the fit. These Ansatze are motivated by the observed rapid temporal 
decay of these functions, modulated by damped oscillations. 



To 


r fitted 


t Eq. jsnj 





0.785418 


0.785398 


0.1 


0.90101 


0.901016 


0.2 


1.04403 


1.04403 


0.3 


1.22641 


1.22641 


0.4 


1.46803 


1.46803 


0.5 


1.80465 


1.80464 


0.6 


2.30773 


2.30771 


0.7 


3.14395 


3.14393 


0.8 


4.81334 


4.81332 


0.9 


9.81476 


9.8159 



TABLE I: Comparison between the characteristic time r extracted by the best fit of Eq. (|84p to the numerical data and the 
corresponding value predicted by Eq. (|86|) for various values of To- The fit has been done within the time interval t £ [20, 60] for 
the correlation function C+(t). We found good agreement between the analytical expression (|86|l and the fit of the numerical 
data also for the decay of the response function. The major source of uncertainty in this fit is the systematic error associated 
with the choice of the interval within which the data are fitted by Eqs. (|84|) and (|85[) . In the worst case, this choice affects the 
estimate of r in its third digit. 



Remarkably, the values of r that provide the best fit to the numerical data turn out to coincide with great accuracy 
(see Table HI with the time constant 

_! r dfc de fc (r) 4arctan(yT^T) 

T = / — — lncosA fe l,r = -== , 86 

Jo T dfc ttVT - 1 

calculated in Refs. [2l|, [23|, which characterizes the exponential long-time decay of (af(t)) oc exp(— t/r) and conse- 
quently of the equal-time correlation (af(t)a x (t)) oc exp(— 2t/r) at large spatial separations \i — j\ At. Although 
we are not aware of any analytical proof of this fact, we conjecture that the time constant r of the long-time expo- 
nential decay of (af(t + to)af(to)) oc exp(— t/r) in the stationary regime to — > oo is determined by Eq. beyond 
the numerical coincidence shown in Table U [7l| . Indeed, the analytic expression in Eq. (|86|) finds further support 
from the fact that it correctly reproduces the characteristic time r oq [72| which controls the long-time decay of the 
equilibrium correlation function (erf (t)df (0)) at temperature T when, according to the formal mapping in Eq. (|43l) . 
one substitutes cosA^ in Eq. (|86[) with tanh(e/ £ /(2T)). In addition, note that the non-equilibrium coherence time r 
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given by Eq. (|86|) decreases upon increasing |1 — Tq\, i.e., the energy injected into the system and r ~ 1 1 — To | 1 
for r — > 1(= r). This fact is consistent with the equilibrium behavior r cq ~ 1/T — > oo for T — > 0: finally, with 
this assumption, the relation between the long-time decay of (af(t)af(s)) for t = s (studied in Refs. [HI, and for 
t, s — > oo with fixed |t — s| (present study) is the same as the one generally expected for the dynamic correlations 
($(r, i)$(0, s)) of the order parameter $(r, t) at point r and time t after a quench to a quantum critical point of a 
one-dimensional model with a linear dispersion relation tk <x |fc| (i.e., a 1+1-dimensional CFT) [3(1]. Analogously to 
Ref. [36| . also in the present case we expect the non-linearity and the upper bound of the dispersion relation to be 
responsible for the oscillating corrections in Eqs. ([53)1 and ([8"5"j) . 




FIG. 10: Decay of the correlation function C% (black line) and of the linear response function R x (green line) of the order 
parameter, as a function of time t for To = 0.3. Upper inset: zoom into the long-i decay that demonstrates the exponential 
relaxation with the characteristic time r defined in Eq. (|86[) (dashed yellow line). Lower inset: (e 1 ^ C% /Ac — l)/oc (black) 
and (e*' r R x /Ar — l)/an (green) vs. t. The yellow dashed line represents the t -1 / 2 envelope of the damped oscillations, in 
agreement with Eqs. (|84)l and (|85|) . 



Figure HU1 presents C+(t) (black) and R x (t) (green) as functions of time t, for To = 0.3. The upper inset shows 
a zoom of the main plot for large times and compares it with the leading exponential decay (yellow dashed curve) 
with the rate r _1 given by Eq. ([86]) . The lower inset shows the correction to the leading decay, (e t / T C 3 j r /Ac — l)/ac 
(black) and (e t ^ T R x /An — I) /an (green), estimated in Eqs. ([53|) and (|55"j) and the yellow dashed line represents the 
envelope t~ x / 2 . Although several fitting parameters are involved in Eqs. and (|85p . we tested these expressions for 
several values of Tq £ {0, 0.1, 0.2, . . . , 0.9} and they turned out to be always remarkably accurate. The exponential 
decay sets in very soon, usually already after the first oscillation of and R x . The next-to leading order oscillatory 
correction involving a power-law decay is more difficult to determine, but still we found a good agreement quite soon, 
with ac — 0.21 and an ~ 0.52 rather independent of the value of Tq. While the parameters A^c depend on Tq, their 
ratio 



Ac/A R = 1.210(5) (87) 

does not within our numerical accuracy and the range considered. In fact, for the fixed initial time to — 10 that we 
have primarily considered in our computations, only the ratio corresponding to Tq = 0.9 is slightly larger than the 
others. We argue that this is due to non-stationary contributions and indeed, upon increasing to, the corresponding 
value of the ratio decreases towards the value obtained for smaller Tq. This observation agrees with the fact that for 
the correlation functions C± to be stationary after the quench, one has generically to take t.Q 3> r, with the coherence 
time t given in Eq. (|56"|) . As this coherence time diverges for To -> 1, the investigation of shallow quenches with 
To ~ r = 1 necessarily requires studying increasingly large values of tQ. In addition, a large coherence time makes the 
numerical determination of the effective temperature in the uj — > limit more difficult as it requires the integration of 
the correlation function over a longer time interval. This is the reason why we have limited our analysis to sufficiently 
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deep quenches with r < to = 10, i.e., from Eq. (J86J) , To ^$ 0.9. In the left panel of Fig. [TT]we show the exponential 
decay of C+(t) for several initial conditions T and its comparison with the rate in Eq. The exponential decay 

of R x (t) is precisely the same. In the right panel we show the damped oscillations displayed by C^_(t) ci /Ac [where 
the value of Ac is determined by the best fit and r from Eq. (|86p] and with a dashed line the envelope oc t^ 1 / 2 . The 
inset of this panel highlights the phase shift </> among the various curves. Note that — apart from this tiny shift — the 
curves almost overlap, suggesting that the oscillatory corrections in Eqs. ([84]) and (|85[) might be actually independent 
of T and that these observed differences are finite-size or numerical effects. 




FIG. 11: Left panel: Time decay of the correlation function C+(t) (solid lines) for several initial values of To. For each of these 
values the dashed lines correspond to an exponential decay with the rate given by Eq. (|86[) . which clearly capture the long-time 
behavior of C+(t). The amplitude Aq of this exponential decay (see Eq. (|84|) ~) is estimated by the best fit. Right panel: 
Oscillatory correction to the leading exponential decay of the correlation C+(t). Solid lines indicate Q, x (t) = C+^e*^ /Ac [see 
Eq. (|84|) ]. with r determined from Eq. ()86p and Ac as indicated above. The dashed line is an envelope 1 ±0.23/%/! that shows 
the amplitude of the damped oscillations. The inset zooms in the time interval t £ [24, 27] and shows the slight phase shift 
among the different curves. In both panels red, green, blue and yellow lines correspond to To = 0, 0.3, 0.5, and 0.8, respectively. 
For the sake of simplicity we do not show here R x (t), but we found the same characteristic time and a similar oscillatory 
behavior. 



Since the correlation and response functions are computed on a chain of finite length L, they display the L- 
independent behavior described by Eqs. ([84]) and (|85j) only for times t smaller than a certain cross-over time set by L, 
after which finite-size effects dominate. Accordingly, in comparing numerical data with these asymptotic expressions 
we restricted to suitable long times within the " early" regime, being the detailed discussion of finite-size effects beyond 
the scope of our study. 

The discussion in Sees. IIV Al and IIVBI reveals that the qualitative behavior of the various effective temperatures 
for the transverse magnetization can be affected by the spatial structure of the quantity under study, as indeed local 
and global quantities display different features. In order to understand the extent up to which the spatial structure of 
the correlation and response functions influences such effective temperatures, we investigated this issue numerically 
by computing the two-time and two-point functions of the local order parameter: 

CI (r, to+t,t Q ) = \{ [tt+r (to + t), of (t )] J, 

which provide the symmetric correlation function Cf and the linear response function R x {r,t$ +t, to) = 2iCZ(r, to 
t,to), for t > [53l ]. The computation of this quantity is still based on the same determinant equation as Eq. 
where the indices j\ and l\ run from | + 1 to L - r, j'2 and I2 from h + 2 to L — r + 1, j'3 and I3 from 1 to -| — 
and J4 and I4 from 2 to 4 — r + 1. We fixed the length of the chain L = 1000 and the waiting time to = 10 and we 
studied the dependence on t after quenches to the critical point T = 1, starting from an initial condition which is the 
ground state of the Hamiltonian corresponding to To € {0, 0.1, 0.2, . . . , 0.9}, exactly as done before for the case r = 0. 

The results of the numerical calculation of C^.(r = 10, t) and R x (r = 10, t) are presented in Fig. [12] (left and right 
panel, respectively) as functions of the time separation t at which the two observables are measured, see Eq. (|88|) . 
Both correlation and response functions are enhanced for small values of 1 1 — Tq | and, as expected, they vanish in 
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C + (r,t) 




FIG. 12: Time-dependence of the stationary correlation and linear response function of two a x spins (order parameter) separated 
by a distance r = 10 in units of the lattice spacing, after a quench to the critical point F = 1. The correlation function 
C+(r = 10, t) is reported in the left panel, whereas the response function R x (r — 10, i) in the right one. From bottom to 
top, the various curves refer to the different initial conditions To = (red), 0.3 (green), 0.5 (blue) and 0.8 (yellow). Both 
C+(r = 10, i) and R x (r = 10, t) display clear light-cone effects, which are discussed in the main text and summarized in Fig. 1141 
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FIG. 13: Time-dependence of the correlation and response functions C%(r, t) (left panel) and R x (r,t) (right panel), respectively, 
for To = 0.3 and, from top to bottom, various values of r = 5 (red curve), 10 (green) and 20 (blue). The behavior of both 
functions at short times is compatible with a light-cone effect with characteristic time r/v m = r/2 for the present case r = 1 
(see the main text). Before this characteristic time, the correlation function is almost constant whereas the response function 
is negligible. The eventual exponential decay (highlighted by the choice of the logarithmic scale) is clearly independent of r 
and the dashed lines correspond to a decay rate given by Eq. (|86|l . which we conjecture to capture the behavior for r = 0. The 
dependence on r of the correlation function C+ (r, t ~ 0) at small times is compatible with a spatial exponential decay with 
correlation length £ — 2.3 given by Eq. (189p , 



the limit of large time separations. However, differently from the case r = shown in Fig. [TUl these correlations with 
r / display a light-cone effect due to the finite maximum speed of the quasi-particles of the model which move 
ballistically, with a maximal speed v m . This effect manifests itself in the fact that both the (connected) correlation 
and the response functions in Fig. II 21 remain almost constant up to times t ~ r/2. After that, the correlation function 
decays oscillating towards its asymptotic vanishing value, whereas the response function first abruptly takes non- 
negligible values and then decays as well. This feature is highlighted in Fig. [T31 which shows on a logarithmic scale 
the functions C+(r,t) (left panel) and R x (r,t) (right panel) with fixed Tq = 0.3 and, from top to bottom, various 
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values of r — 5 (red), 10 (green), 20 (blue). The behavior can be understood by extending to dynamical quantities in 
the stationary case the qualitative picture of Refs. [U |H, HjJ as depicted and summarized in Fig. Q3] Indeed, 
one expects the response function R x (V, t > 0) to be negligible (vanishing in the scaling limit in analogy with what 
happens for the correlation function [23[) and the correlation function C+(r, t) to be almost constant as long as the 
point located at (r,t + to) does not belong to the forward "light cone" of the point (r = 0,to). Taking into account 
that the forward light cone expands with a speed w m = maxt |defc(r)/dfc| = 2min{r, 1} around the vertex point, this 
is compatible with the horizon effect that appears for T = 1 at t ~ r/v m = r/2. Note that the emergence of a finite 
speed at which correlations and responses propagate is also present at equilibrium [73[ . 



time 




! v mtn ! v m s 1 

I U 1 

FIG. 14: Schematic plot of the various light-cone effects after a quench occurring at t = 0, in the "space-time" diagram. 
Correlations between two points A and B separated by a spatial distance r change because of the quench only if their backward 
light-cones, indicated in the figure by dashed lines originating from the points, include spatial points which were significantly 
correlated in the initial state. This is not the case for A and Bi (assuming the spatial extent £o of the initial correlations to be 
negligible compared to the distance between the cones) but it occurs for A and B2, as there is an overlap in the point S, which 
extends into a wider region at later times [3(|. Indicating by to and s the times of A and B2, respectively, the existence of an 
overlap requires v m (to + s) > r, where v m is the speed with which the light cone extends around a point. For an equal-time 
correlation function such as the one between A and B3 (with s = to), this occurs only if a time larger than the "Fermi time" 
tp — r/(2u m ) [23|] has elapsed since the quench, while equal-time connected correlation functions are approximately constant 
and equal to their equilibrium values at zero temperature for to < tp. Similarly, the symmetric correlation function is expected 
to be constant and the response function negligible if the point B does not belong to the forward light cone (indicated by the 
shaded area) of the point A at which the perturbation is applied. This occurs if t ~ r/v m = 2 x tp, as it can be seen graphically 
from the figure and as discussed in the main text. 

The eventual exponential time-decay of both C+(r,i) and R x (r,t) in Fig.[T2]is actually independent of r and occurs 
with the rate in Eq. (|86| which we conjectured to describe the exponential decay of the correlation and response 
functions with r = [see Eqs. (|84[) . (pj|) and Table H], indicated by the dashed lines in Fig. [T3] Focusing on the 
behavior of the correlation function C+ at short times, it is possible to study the dependence of its plateau values for 
t < r/v m = r/2 on the distance r and in particular how they decay as functions of r. This issue has been addressed 
in Refs. [2l],[23| for (af +r (to)af (to)) [which corresponds to the case t = of the dynamic correlation C+{r,t) studied 
here] and, as expected, the spatial de-corrclation of these plateau values in the present case follows an exponential 
decay with the correlation length £ found in Refs. pll. [23j: 

C' = -£^ lncosA fc (l,r ) = ln(l + -l=) . (89) 
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3. Effective temperatures 



Here we focus on the effective temperatures associated with the stationary dynamics of the order parameter a x both 
in the time and in the frequency domain at small frequencies u> — > 0, starting from the latter. Due to the numerical 
nature of the calculation of C±(t) as a function of time, it is very difficult to access reliably the full frequency 

dependence of C±(uj). Still, the numerical results obtained for cj — > are robust, also against the fact that a finite 
time interval for the integration has necessarily to be taken. In fact, we have studied the small- frequency regime in 
a time interval t G [0,£ max ] for different values of t max and we found that it stabilizes rapidly with increasing i max , 
t mas being an increasing function of t. Small frequencies capture the exponential decay and they are less sensitive 
to other details. For these reasons, their extrapolation in a sufficiently large, but finite, time interval is expected to 
provide reliable results. Even if a complete characterization in the full frequency spectrum of the correlation functions 
and the associated effective temperatures is definitely an interesting issue, the region of small w is for us the most 
interesting because we aim at comparing P x ^{^ — > 0) with /3* ff obtained from Eq. (|20[) in the long-time limit. The 
zero-frequency limit of Eq. (|19p yields: 

f™dttR x (t) 

that we take as our working definition of (3 X ff(t<J — > 0). The corresponding temperature T£ s (lj = 0) = 1//3® s (lj = 0) 
is indicated by circles in Fig. [15] as a function of To < 1. Alternatively, as discussed after Eq. (|21|). we can enforce 
a constant value /3* s in the (generalized) FDT (J2UJ) in the time domain. This operation allows one to interpret the 
r.h.s. of that equation as a series of time derivatives of C+(i) and, after the substitution into this series of the long-time 
behaviors in Eqs. ([84]) and (|85j) . it yields 

HA R /(2A C ) = tan(fy? c V2r) (91) 

for t — > oo, in which one neglects the oscillatory corrections in Eqs. (|84p and (|5S"|) . [Note that the oscillatory terms 
could not be discarded in the analogous analysis presented for M after Eq. ([74]) because they were actually providing 
the leading contributions.] Note that h has been reinstated in Eq. (|9T|) for completeness, but we shall set again h = 1 
in what follows. As we pointed out in Eq. (|87| . Ar/Ac — and therefore the l.h.s. of Eq. (|9Tj) — does not depend 
significantly on To within the range of values of Tq investigated here; accordingly, (3 X S on the r.h.s. of Eq. (f9"Tj) can 
be expressed in terms of r as /3 X S ~ 0.78t and inherits its dependence on Tq given by Eq. (|8<3|) . The corresponding 
temperature T x s = 1//3 X S is referred to as il T x s t ^> 1" in Fig. [15j where it is plotted as a function of r and the 
corresponding data points are denoted by diamonds. In this case we find that the effective temperatures determined in 
frequency (circles) and time domain (diamonds) are almost indistinguishable, especially for To ~ 1. For H(3 x s /2t <C 1 
in Eq. flSTJ one recovers the classical limit {23) /3* ff ~ -R x (t)/[dC x (t)/dt] ~ tA r /A c ~ 0.83r. The value of the 
corresponding temperature, referred to as ll T x s t ^> 1 FDT class." is indicated by squares in Fig. [151 All the three 
determinations of T x s discussed so far are compared in Fig. [15] with the effective temperature T^ s (dashed line) 
determined on the basis of Eq. (|4"5"]) [30| . In the same figure we also show (triangles) the zero- frequency effective 
temperature obtained from Eq. ([90]) on the basis of the two-point functions C+(r,t) and R x {r,f) with r = 10. 

In order to facilitate the quantitative comparison between the various effective temperatures T" s (which we label 
generically by a) reported in Fig. [T5] its right panel shows their relative difference AT" = (X^ — T^ s )/T^ as a 
function of To, where the effective temperature Tj^ (dashed lines in both panels of the figure) obtained from the 
energy balance is taken as the reference. As we discuss in App. [D] 



rpE „ 2 y/6 (logT + 41og2-3) 



for T — > oo, i.e., in the limit of shallow quenches ro —> T = 1 (see Eq. ([54])). Accordingly, taking into account 
Eq. (J7U), the ratio AT M = (T^ - T^ s )/T^ s tends to -1 for T -> 1, with the logarithmic approach clearly displayed 
by the dash-dotted curve on the right panel of Fig. [15] Both panels of this figure demonstrate that, in general, the 
temperatures defined from the FDR associated with the stationary dynamics of the various quantities do not coincide 
with the "static" temperature Tj^ defined on the basis of the energy of the system. However, in the neighborhood of 
Tq ~ 0.4, the relative discrepancies ATjfg are less pronounced (see right panel), which can be heuristically traced back 
to the fact that the temperature T^ s of the various modes (see Eq. ([73]) and the right panel of Fig. [7]) is practically 
constant and equal to within a rather wide interval of momenta k. Upon moving away from Tq ~ 0.4, such an 

interval shrinks, as it is clearly shown by the various curves on the right panel of Fig. [7] which correspond to different 
values of Tq. On the other hand, differently from the cases of M and a z , the frequency-dependent dynamic effective 
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FIG. 15: Left panel: various parameters TJj (symbols and solid lines, see the key) as a function of To, compared with the 
effective temperature T e g defined from the energy [dashed line, see Eq. (|45[) ]. The solid lines, from bottom to top, indicate 
the values determined on the basis of the classical limit of the FDR in the time domain (squares), of the limit ui — > of the 
frequency-domain FDR (circles), of Eq. (|91[) (diamonds), and of the limit u — > of the frequency-domain FDR but for spins 
separated by a distance r = 10 (triangles). The orange dashed-dotted line shows the effective temperature T^g (lj — > 0) obtained 
for the total transverse magnetization in the limit of w — > 0. Right panel: Relative differences AT a = (T" H — Tjj)/Tj^ between 
the various effective temperatures reported on the left panel and T eS taken as a reference, as functions of Fo. Colors and 
symbols are the same as those of the left panel. 

temperature for the order parameter takes a non- vanishing value for to = (circles in Fig. I15[) , which can be recovered 
with great accuracy from the long-time limit of the (classical and especially quantum) FDR in the time domain 
(squares and diamonds in Fig. IT5|) . Whereas these values are almost indistinguishable for To ~ 1, slight discrepancies 
emerge away from the critical point. However, these might be due to the numerical accuracy of the calculation. More 
consistent discrepancies, instead, emerge in comparison with the case r = 10 (and with the temperature T^), but 
they become anyhow negligible for 0.5 < Tq < 2, i.e., in a rather wide neighborhood of the critical point. 

V. SUMMARY AND CONCLUSIONS 

In this Section we first summarize our findings and we then discuss the perspectives of this work, which was partly 
anticipated by the brief account in Ref. [f|. 

Aiming at understanding whether and how thermalization to the canonical ensemble described by a Gibbs distri- 
bution may arise after a quench of an isolated quantum system, we focused on the fluctuation-dissipation relations 
(FDRs) between a set of dynamic quantities. By comparing these FDRs to the canonical FDT holding in equilibrium, 
we extracted parameters, actually functions, that in any (even partial) equilibrium situation should be equal for all 
choices of observables and, moreover, constant. With a certain abuse of terminology we called these functions effective 
temperatures and we examined whether they satisfy, and under which circumstances they do, the constancy property 
just mentioned. 

For illustration purposes we pursued this approach by investigating the dynamics of the quantum Ising chain in 
a transverse field T (see Sec. IIII[) . which is initially prepared in the ground state of the Hamiltonian H(T) with 
r = Tq and then quenched to the critical point r = 1. We computed correlation and (linear) response functions of 
the local transverse magnetization <r z , the global magnetization M z = (J^. of )/L and the local order parameter a x , 
extracting the associated effective temperatures. We chose to work with this model knowing that sound evidence has 
accumulated over the last years for the fact that its stationary properties, as well as those of more gen eral integrable 
systems, can be consistently described by the so-called generalized Gibbs ensemble (GGE) (l7l. [iH, . This implies 
that the effective temperatures defined from the canonical FDRs are not expected a priori to take a common value. 
Nonetheless, they may anyhow signal regimes of partial equilibration and thus capture important features of the 
dynamics of the system, as it happens in a number of classical non-equilibrium cases in which the Gibbs stationary 
ensemble cannot be reached because of the slow relaxation [l4| • 

As a first observation let us note that the FDRs can be used to identify effective temperatures independently of 
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the functional form of the correlation and (linear) response functions involved in the FDRs. This is of great practical 
advantage since it allows us to investigate also quantities with power-law (or any other kind of) decays, as opposed 
to the strategies that base the comparison between the quenched (non-equilibrium) and the thermal (equilibrium) 
dynamics on the analysis of the characteristic time- or length-scales primarily associated with exponential decays. 
Moreover, it is well-known from the analysis of classical [l2| (and to a certain extent also quantum [43|, [5l|) glassy and 
coarsening dissipative models that the functional form of the decay of the correlation functions can be highly non- 
trivial (non-stationary, non-exponential) and yet, in some circumstances, an effective thermal behavior can develop 
asymptotically in certain time regimes. It is therefore highly desirable not to confuse the qualitative features of time 
decays with the possible emergence of a thermal behavior. The investigation of the system's properties in terms of 
FDRs provides the possibility to analyze separately its behavior within different time (or frequency) regimes. 

In order to illustrate this general approach, we calculated the (self) FDR for three observables of the Ising model 
that are local (a* ' z ) or non-local (M) in space and local (of, M) or non-local (of ) in the quasi-particles. By inspection 
of the FDR and of the frequency-dependent T c z ff parameter that we extract from it, we concluded that the dynamics 
of of is not compatible with Gibbs thermal equilibrium at any effective temperature, in spite of the fact that thermal- 
like behavior — compatible with a Gibbs distribution at an effective temperature Tj^ set by the initial energy - 
is observed for the stationary expectation value (o~ z )q = (<7 z ) t=T -b . Note that this instance clearly demonstrates 
the importance of the alternative approach to the issue of thermalization we are currently proposing beyond the 
specific model investigated here: indeed, the analysis of specific time- independent quantities (in this case, (ct z )q) 
would suggest a picture of thermalization which can be easily disproved by focussing on dynamical quantities (in this 
case, the two-time response and correlation functions <r z ). 

The FDR in the frequency-domain yields for M a finite T^g(To) in the limit ui — > (see Fig. [7J right panel and the 
dot-dashed line in Fig. [8]) , which is actually connected to the fact that low-energy modes have a finite mode-dependent 
effective temperature T e ^ [see Eq. (|35))] for k — > 0, in relation to the GGE [see Eq. ([75]) ]. This feature is peculiar 
of quenches to (and from) the critical point and it does not carry over to the other cases. However, contrary to the 
heuristic expectation, we did not find a way to recover the value T^(Tq) directly from the long-time limit of the FDR 
in the time domain. In fact, the effective temperature T* s of the long-time dynamics of C±(t) (and derivatives) is 
somehow controlled by the behavior at high frequencies, as discussed in Sec. IIVBI This is a natural consequence of 
the presence of oscillating terms at the leading order in C± (t) with a frequency equal to the frequency w max above 
which the Fourier transform of C± (t) vanishes identically. This structure highlights that not only the slow power-law 

decay of C±' z (t) is important, but also the oscillations are a constitutive ingredient. This structure, as well as the 
presence of the maximum frequency w max , is the combined result of the presence of the lattice cut-off that bounds the 
dispersion relation and of the particular quadratic form of the observable with respect to the elementary excitations 
of the model. 

The dynamics of the operator of shows, instead, a very different behavior. Both C+(t) and R x (t) decay expo- 
nentially (see Figs. [TU1 ITT1 [TBI and the discussion further below) with a characteristic time given by Eq. (|55|) . which 
decreases upon increasing the energy injected into the system. Differently from the correlations of a z , the time decay 
of C+(t) and R x (t), for quenches at the critical point T = 1, does not display oscillations at the leading order for 
large times. The oscillatory part, which intervenes only in the leading corrections, is therefore less relevant and one 
would expect that the long-time limit of the corresponding T* ff could be regularly recovered from a small-frequency 
expansion, as it is actually the case (compare circles and diamonds in Fig. [Tip]) . and 23r have the same qualitative 
dependence on To but they still present some differences (and they also differ from Tj^), thus excluding a single 
temperature effective Gibbs description. Finally, we have shown that, by choosing appropriate observables, one can 
recover from the study of FDRs the effective temperatures T e fc ff associated with the modes, which in turn characterize 
the GGE [see Eq. ([75]l]. 

While the primary focus of this work was on the application of fluctuation-dissipation relations for addressing the 
issue of the possible thermalization of isolated quantum systems long after a quench, our study of two-time quantities 
(connected correlation and response functions) lead to the following observations as far as their behavior is concerned: 
(i) For a quench to the critical point T = 1, the dynamical two-time quantities after the quench are invariant under 
To <H- r^ 1 in the stationary regime. Such an invariancc is broken when at least one of the two times approaches the 
moment of the quench. This statement is supported by analytical calculations in the case of the local and global 
transverse magnetizations of (see App. [A} and M z , whereas it is based on numerical observations in the case of 
of. Moreover, the expectation value of the energy [see Eq. pi))) ] and of the transverse magnetization [see Eq. (|47|)] 
for r = 1 possess the same invariance as functions of Tq. We conjecture that this should be not only a generic 
property of the dynamics in the stationary state but it should also be largely independent of the specific quantity 
under study, being somehow connected to the duality of the equilibrium model. To us, a deeper understanding of 
the nature and the limits of such a symmetry — which is broken by non-stationary contributions — is definitely 
an interesting issue, (ii) In the stationary state, the correlation and response functions of the order parameter of 
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decay exponentially for large time separations, as a function of it [sec Eqs. (|84[) and (|85j) ]. The decoherence time r 
of this decay [see Eq. ([55)1 ] is (numerically) the same as the one which characterizes the relaxation of the average 
of the order parameter (of (t)) oc exp(— t/r) after the quench from the ferromagnetic phase and of the equal-time 
correlation (of (t)crj(t)) oc exp(— 2t/r) for \i—j\ 3> 4t, calculated in Refs. [2l|,[23|. Even though we only have numerical 
evidence for this fact (for T = 1), we conjecture that this should be true beyond this case, e.g., for quenches towards 
non-critical states (see Sec. II V CI for additional arguments). An analytical proof of these features [and particularly 
of Eqs. and (|85p] should be possible within the approaches discussed in Refs. [U, [2^, [68| . (iii) The two-point 
dynamical correlation functions of the order parameter in the stationary regime exhibit — similarly to the relaxation 
dynamics of the two-points and equal-time correlations [2l], [2^, |45|, [46| — light-cone effects which are compatible with 
the expectations based on the results for correlation functions of CFTs (i.e., in the scaling limit) (36J. This picture 
turns out to extend also to the response function, which is negligible up to times t ~ r/v max =r/2 (see Fig. 1141 for a 
comprehensive schematic summary of these effects). 

In conclusion, concerning the issue of thermalization in quantum quenches which primarily motivated our study, 
we emphasize that a bona fide thermal behavior should always be accompanied by the validity of suitable FDRs, 
also in the context of isolated systems. In this respect it would be desirable to compute FDRs in the cases in which 
thermalization to a Gibbs ensemble is eventually expected, such as non-integrable systems. However, due to the 
general lack of analytic solutions even for the simplest models, this would require the use of numerical methods, 
for which it is still rather difficult to access the long-time stationary regime. In spite of this difficulty, it is worth 
mentioning that the approach we propose here for probing the possible thermalization of a quantum system involves 
only basic quantities — correlation and response functions — which, at least in principle, can be naturally accessed 
both in numerical and experimental investigations. Drawing an analogy with the case of non-equilibrium classical 
systems, the effective temperatures defined above might still provide a useful description of some coarser aspects of 
the physics of the system or of its thermodynamic, such as, e.g., its energy exchange with or response to a device used 
as a thermometer. Clearly, this question is of interest especially when there is no eventual thermalization to a Gibbs 
state, e.g., in the case of the Ising model considered in the present work or more generally for integrable systems. 
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Appendix A: Transverse magnetization 

This Appendix provides the details of some analytical calculations concerning the transverse magnetization whose 
results have been presented in the main text. As explained in Sec. IIII Al we consider a chain of finite size L and 
then we take the thermodynamic limit, which formally amounts to the substitution 1/L y^. — > dk/(2n) j3l| . In 
the following we focus on the (connected) symmetrized and antisymmetrized autocorrelation functions C+ and C_, 
respectively: 

C z + (t + t , t ) = i (Vol R? (f + t ) , of (to)} |Vo) - (Vol of (t + to) |Vo) (Vol of (t ) IVo) (Al) 

and 

C z _(t + t ,t ) = 5(^b|[of(i + to),of(t )]|^)) J (A2) 

where the expectation value is taken over the ground state |Vo) = |0) r of the Hamiltonian H(Tq) before the quench. 
We consider the connected correlations since in general the expectation value (Vol of (t) IVo) does not vanish. The 
Kubo formula [see Eq. ([12]) ] allows one to express the linear response function R z in terms of C_, whereas Eq. (|T3]) 
connects C+ and C_ to the real and imaginary part, respectively, of the correlation function C z (t+to, to) = (Vol of (t+ 
to)of (to) |Vo)- As outlined in Sec. IIII Bl this expectation value is calculated by first expressing of (t) in terms of Ci(t) 
[see Eq. (|28|) ] and then Cj(t) in terms of the quasi- particles {^1°} of H(T ) via the matrix elements defined in 
Eq. (|37[) . Introducing the simplified notation Ufc(t) = v^' r °(t) and Mfc(t) = w^ ,r °(t), the symmetric and antisymmetric 
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autocorrelations for the transverse magnetization are given by the real and the imaginary part (see Eqs. (|13|) and 
O) of 



C 



(t + ta,t ) = JoYj [vk(t + to)v* k (t )ui(t + to)u^{t ) + v k {t + t )vt(to)ui{t + t )u* k (t ) . (A3) 

k,l 



Note that, according to Eqs. ([35} and Q5JJ) , 

u_fe = Ufc whereas v- k 



-Vk- 



(A4) 



This implies that the second term in the sums is odd with respect to k and therefore it does not contribute to the total 
(symmetric) sum. Then, expressing the first term via the angles of the Bogoliubov rotation introduced in Eqs. (|38j) 
and (|37[) . one eventually finds (e k = £fc(r), see Eq. 



C z + {t + t , t ) V] [cos(e fe t) (sin 2 9 k cos 2 4 + cos 2 6 k sin 2 4) - 2cos(e fc (i + 2t )) sin 6 k cos 4 cos 4 sin 4] 



[cos(e;t) (cos 2 0; cos 2 5; + sin 2 9i sin 2 5;) + 2 cos(ef (£ + 2io)) sin^ cos Si cos 9i sin<5j] 
~J2 ^2 [ sm ( £ fcO (sin 2 9 k cos 2 S k - cos 2 4 sin 2 4)] ^ [sin(ezi) (cos 2 9i cos 2 <5; - sin 2 9i sin 2 <5;)] 



(A5) 



In this expression (and in the analogous one for Cf, which can be obtained from the imaginary part of Eq. (|A3|) ) one 
can easily recognize a stationary part C+(t) = limt^oo C+(t + to, to) which depends only on t and a non-stationary 
contribution Cf ,ns (£ + to, to) = C+(i + ^o, to) — C+ft ) which depend on both times t and to- We anticipate here that 
the long-time limit to — > oo of the terms in Eq. (|A5|) which involve cos(efe(t + 2to)) vanishes, as it is clear both from a 
numerical study of the corresponding sums and from the analysis of the long-time behavior of similar expressions in 
the thermodynamic limit, done further below in App. IA 31 (The same applies to the identical terms which appear in 
the analogous equation for Cf:(t + to,to), which, however, we do not report here.) Accordingly, the stationary part of 
is given by 



C z + {t) =j2^Z [ C0S M) ( sin2 °k cos2 5 k + cos 2 9, 

k 

4 

~ Z 2 



fc sin 2 4 



cos 2 9\ cos 2 5i + sin 2 9\ sin 2 5i 



Y [cos(e/t) 
i 

^ ^sin(efet)^sin 2 9 k cos 2 4 — cos 2 9 k sin 2 4^ J ^ [sin(e;f) ^ cos 2 9\ cos 2 Si — sin 2 9i sin 2 Si 



(A6) 



For later convenience we report here also the expression of the non-stationary contribution C+ ns (t + to, to), which 
takes the form: 



C+ n8 (t + t , to) X/ [ cos ( e fct)(sin 2 9 k - cos 2 4) (cos 2 4 - sin 2 4)2cos(e;(t + 2t )) sin6> ; cos (5; cos6> ; sine 

l,k 

— ^ 2^ 2 cos(e/ (t + 2t )) sin 9\ cos 5; cos 9\ sin <5; 
f 

= -\Y^ [cos(e fc t)cos(24)cosA fe ]t¥ + (t,to)-^ 2 (t,t ) 



(A7) 



where we introduced A k = 24, 



W+(t,f ) = Y^cos(e;(t + 2t ))sin(26'0sinA i 



and 



sinA fe (r,r ) 



(Tq — r) sin k 



+ r 2 - 2r cos ky/i + r 2 - 2r cos k 



(A8) 



(A9) 
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As expected from their very definitions, W+{t,to) and C+ ni! decay to zero for to —> oo. The stationary part Cl(i) 
limto^oo Cf. (i + i 0l *o) 01 C- + *o> *o)j instead, is given by 



C* (t) [sin(e**)(- sin 2 fc + sin 2 tf fc )] ^ { cos(e^) [cos 2 0, - cos 2 (20, ) sin 2 5,] } 

fc i 

+ £2 ^2 { cos ( efei ) sin2 ^ + cos2 ( 2 ^) sin2 S k] } [sin(e;i)(- cosfy 2 + sin 5?) 



(A10) 



Due to the fact that the terms in the sums of Eqs. (|A6[) and (|A10[) are all even functions of k, in what follows we will 
replace ^ fe with 2^ /c>0 . Note that the expressions reported above apply to an arbitrary quench. 



1. Quenches to the critical point T = 1 



The primary interest of the present study is in the case of a quench to the critical point T = 1, for which we can 
further specialize the equations reported above. In particular, the r.h.s. of Eq. (|32j) is always positive, and therefore 
this relation for k > can be inverted with < 29 k < tt/2: 



cos(2#fc) 



1 



V'l + tan 2 (20 fc ) 



sin(fc/2). 



(All) 



In view of this equation and of Eq. ([321) for cos Afc one can conveniently express all the trigonometric functions of 
Sk and 9k in Eqs. (|A6|) and (|A10[) in terms of cos(25fc) and cos(20fc). Indeed, C±(i) in the stationary state depend 
on Sk — Afc/ 2 and 9k via sin 2 and cos 2 and therefore they can unambiguously be expressed in terms of the cos of 
the corresponding double angles, i.e., A& and 29k, and in particular in terms of cosAfc(r,r ), which encodes the 
dependence on the initial state. However, this is not the case for the non-stationary contribution, which indeed 
requires also the introduction of the sin of these angles, as in Eqs. (|A7|) and (|A8|) . and in particular of sin Afc(r,r ) 
(see Eq. (|A9[) ). As a consequence, the stationary and non-stationary parts of C± at the critical point r = 1 display 
different behaviors under the mapping To — > Tq 1 of the initial condition. In particular, as we discuss after Eq. (|54[) in 
Sec. HVAl A k {T = 1, T ) = -A fc (T = 1, I^ 1 ) and therefore cos A fc (r = 1, T ) is invariant under this mapping, whereas 
sin A^. (r = l,r ) changes sign, as it can be easily verified also by a direct analysis of Eq. (|A9[) . Accordingly, the 
stationary parts C±(£) of the dynamics in Eqs. (|A6|) and (|A10I) for T = 1 arc invariant for To ^ T^ 1 — a fact which 
will be apparent from the expressions discussed further below — while W(t, to) changes sign, making C+ ns (t + to, to) 
not invariant. According to the general strategy outline in Sec. Ill B I we will focus below only on the stationary parts 
of the correlation and response function, given by Eq. (|A6j) and related to Eq. (|A10j) . respectively. In particular, for 
quenches ending at the critical point T = 1, one has e& = efc(r = 1) = 4| sin(fc/2)| [see Eq. (|34|)]. 



cosA fc (r = i,r c 



[see Eq. (|4"Tj) ]. where A fc = 2Sk and 



2(r + 1) 
e fc (r ) 



sin(fc/2)| 



T|sin(fc/2)| 



\/l + (T-l)sin 2 (fc/2) 



( i±r,, 



(A12) 



(A13) 



[Note that T and therefore cos Afc (r = 1, To) are invariant under the transformation To n- T 1 .] Accordingly, Eq. (|A6 
can be simplified as 



Cl(t) 



dk 



cos(e k t) 



Uo 



C dk 

/ — cos(efei) cos Afc sin(fc/2) 

/O T 



— sin(efci) sin(fc/2) 



f dk . , . 

/ — sm(e fc i) cos Afc 

/o I" 



(A14) 



while the stationary response takes the form [see Eqs. (|12[) and (|A10[) ]: 



R z (t) =4 0(i) 



f* dk , . . f n dl , . 

/ — sm(efct)cosAfc / — cos(e;t) 

Jo 77 Jo n 



f v dk f w dl 

/ — sin(efci) sin(fc/2) / — cos(e;i) sin(^/2) cos A/ 
/o 7T J IT 



(A15) 
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FIG. 16: Dependence of the functions E(t) and F(t) (see Eqs. (|A18[) and (|B7|) . respectively) on the time r for various values 
of Tq. The blue, red, yellow and green solid lines refer respectively to To = 0, 0.3, 0.5, and 0.8. The curves corresponding 
to To > 1 can be obtained from those with To < 1 taking into account that they are both functions of T (see Eq. (|A13j0 
and therefore invariant for To — > Tq . The dashed lines show the corresponding leading asymptotic behaviors for large r (see 
Eqs. (|A36|) and (|B24|) . respectively), which are actually attained rather early in the evolution. The thin dashed lines highlight 
the asymptotic leading algebraic time decay of the envelope ~ 2/\/2tyt. 



By direct inspection of the previous equations, one realizes that and R z can be conveniently expressed in terms 
of the following integrals 



/ — cos{e k t) = J (4t), 
Jo k 



77 dfc 



sin(efct) sin(fc/2) 



* dfc 



sin(e fe i) cos A* 



_l_d 
~l~dt 



/ — cos(ekt) cos Afe sin(fc/2) 
Jo K 



~E(4t)=E'(4t), 



(A16) 
(A17) 
(A18) 
(A19) 



(see, e.g., 10.9.1 and 10.6.2 in Ref. [6J])j where J a is the Bessel function of the first kind and order a, as 

C**(t) = J 2 (4i)- [E'(4t)} 2 + J*(4t)-E 2 (4t), and (A20) 
R z {t) = 46{t)[E(4t)J (4t) - J 1 (4t)E'{4t)}. (A21) 

The function E introduced in Eq. (|A18[) can be explicitly written as 



E(t) 



dk . , . VT sin(fc/2) 

— sin(rsin(fc/2))— v ' ' . 

n ^Jl + {T -l)sm 2 (k/2) 



(A22) 



where T is given in Eq. (|A13|) . The left panel of Fig. [T5] shows E(t) as a function of the time r for various values of 
To. The leading asymptotic behavior of the function for r ^ 1 — which is discussed further below in App. IA 31 — is 
indicated by dashed lines. The right panel of the figure, instead, presents for comparison the behavior of the function 
F which will be introduced in App. [B] [see Eq. (|B7|) ]. Later on in the analysis of the effective temperatures, we will 
need E'(0), which can be readily expressed as 



E'(0) 



2VT E e (l - T) - K e (l - T) 



T 



(A23) 
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where KJk), Ejk) are the complete elliptic integrals of the first and second kind, respectively (see, e.g., chapter 19 
in Ref. [64j|). For later convenience we report here their definitions: 

KJk) = f 1 , d£ -=^=, (A24) 

1 ; Jo VT^Vl-ke 2 V ; 



(A25) 



Compared to the standard notation for these functions (see, e.g., Ref. |64j), we have added a subscript e to the 
corresponding symbols in order to avoid confusion between E and E e . In addition, adopting the same convention as 
Wolfram Mathematica, we indicate by k the square of the moduli of the elliptic functions, which are usually defined 
in the mathematical literature as in Eqs. (|A24[) and (|A25|) with k n- k 2 ; see, e.g., chapter 19 in Ref. fo4| . 

2. Correlation and response functions for critical quenches with initial conditions To = 0, 1, oo 

For To = or To = oo and T = 1, the parameter T in Eq. (|A12[) takes the value T = 1 [see Eq. (|A13p ] and therefore 
cos Afe = | sin(fc/2)|. Accordingly, the function E [see Eq. (|A18[) ] which appears in Eqs. (|A20|) and (|A21[) for C± and 
R z , respectively, can be expressed as a Bessel function (compare with Eq. (|A17|) and see, e.g., 10.6.1 in Ref. |64|): 

E{t) = J 1 {t) and therefore E'(t) = J[(t) = [J (r) - J 2 (r)]/2. (A26) 

Accordingly, 

C%(t) = J 2 (At) - I[J (4t) - J 2 (4t)] 2 , (A27) 

R z (t) = 2 0(f)Ji(4t)[J o (4t) + J 2 (4t)]. (A28) 

The case of the critical quench with To = 1(= T), i.e., T — > oo (sec Eq. (|A13[1 ) actually corresponds to the equilib- 
rium situation at T = 0, for which one expects the correlation and response functions to be stationary immediately 
after the "quench" at t — because, effectively, no quench takes place. The corresponding expressions for and R z 
are readily derived from Eqs. (|A20j) and (|A21j) . by taking into account that in this case cosAfe(r = 1,Tq = 1) = 1 
[see Eq. (|A12[) ] and therefore 

E(t)=H (t), with E'(t)=H-i(t), (A29) 

where H a are the so-called Struve functions (see, e.g., 11.5.1 and 11.4.27 in Ref. [iH). Accordingly, the correlation 
and response functions are 

C% (t ) = J 2 (At) - Hi l (At) + J\ {At ) - Hi (At) , and ( A30) 

R z (t) = A6(t)[H (At)J Q (At) - ,h(At)H^(At)]. (A31) 

The expressions (|A30[) and (|A31[) agree with those for the symmetric correlation derived in Ref . 14411 and for the 
response function after quenches originating from the fully polarized state To = oo derived in Ref. [16( ■ 

3. Asymptotic expansions for large times 

For quenches at the critical point T = 1 and generic value of To, we are interested in the asymptotic expansion 
for t — > oo of the correlation function and of the response function R z in Eqs. (|A20|) and (|A21|) . respectively. 
These expressions involve the function E(t) introduced in Eq. (IA18I) and the Bessel functions Jo.i- The asymptotic 
expansion of the latter is well-known (see, e.g., 10.17.3 in Ref. [641 ]) and therefore we need only to determine the 
asymptotic behavior of E(t), defined in Eqs. (|A18P and (|A22[) . Introducing the integration variable y = 1 — sin(fc/2), 



the latter expression becomes (the procedure below is analogous to the one adopted in Ref. 7J], in particular see 
App. B.2 therein): 

E(r) = I' dy M(l- g )r) x 2VT (1 - y) 

Jo VV * V(2-2/)[l + (T-l)(l-y) 2 ] 
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The asymptotic behavior of this integral can be inferred from the one of 



dy— f(y) 



(A33) 



where f(y) is a generic function assumed to have a regular expansion for y ~ 0, i.e., f(y) = X)^Lo fnV n i such that 
e( T ; /) = *£n=0 /n(-*d/dr) n e(r; 1), with e(r; 1) = v^tt/t x [C{y/2rpK) + iS( y/^r/n)]. Here C and 5 are the so- 
called Fresnel integrals: see, e.g., 7.2.7 and 7.2.8 in Ref. [iH for their definitions and 7.5.3, 7.5.4, 7.12.2 and 7.12.3 
for the corresponding asymptotic expansions. The asymptotic behavior of Eq. (|A33[) for large r can be calculated on 
the basis of these expansions and turns out to be: 



e(r;/) =y/ne 



i-rr/4 



/(0) if'(0) 3/"(0) 
t 1 / 2 2r 3 / 2 8r 5 / 2 



+ 0(t- 7 / 2 ) 



1 



(A34) 



- -2^)^/(1) + -/'(I) + 0(r- 3 )- 



Note that, alternatively, the first line on the r.h.s. of this equation could have been inferred from the stationary-phase 
approximation [67| . One can take advantage of the previous asymptotic expansion in order to calculate the one of the 
integral 



p i(l-y)r 

dy — /(y) = e-e*(r,/) 



/(o) 



sin(r + 7r/4) /'(0) cos(r + 7f/4) 



+ iy/Tr 



-/(o) 



r i/2 2 r 3 / 2 

cos(r + 7r/4) /(l) 1 /' (0)sin(r + 7r/4) 



(A35) 



-1/2 



r 3/2 



+ 0(T" 2 ), 



the imaginary part of which has the same structure as E in Eq. (|A32[) . Taking into account the specific expression of 
/ in Eq. (|A32[) one eventually finds 



E(r) = 



1 



2 7T. 

. "ttt cos(r + -) 



4T - 1 8 • / 7T . . 5/2 . 

— — - sin (t + - + O r- 5 / 2 , 

32V2T r 3 / 2 4 ; 



(A36) 



where the term cx 1/r in Eq. (|A35[) — which emerges upon resumming the expansion of f(u) — drops out because 
/(l) = in the case of a quench with T _1 7^ (which excludes the equilibrium case at zero temperature Tq = 1(= T) 
discussed in Sec. IA 21 and further below, corresponding to T _1 = 0). By using Eqs. (|A36[) and the asymptotic 
expansions of Jq.i (see, e.g., 10.17.3 in Ref. [(34j|) in Eqs. (|A20[) and (|A21|) . we finally find: 



C z + {t) = -—5 cos(8t) + 0(t~ J ) 



R z (t) 



8nt 
1 



Ant 2 



T _1 - sin(8i) + 0(t 



-3\ 



(A37) 
(A38) 



for t 3> 1. In addition, in the special case T = 1 — see Sec. IA 21 — E(t) = Ji(r) and indeed one recovers from 
Eq. (|A36|) the well-known expansion of the Bessel function Ji (see, e.g., 10.17.3 in Ref. (64j). 

In the equilibrium case at zero temperature T" 1 ~ 0, f(l) = 2/ir and Eq. (|A36|) with the additional term 2/(jrr) 
stemming from f(l) 7^ reproduces the known large-T expansion of the Struve function Hq(t) (see 11.6 in Ref. [64j), 
which E reduces to in this case [see Eq. (|A29|i ] . Accordingly, one finds 



cm 
R z (t) 



1 / . K N 

— = cos(4i + — ) 

\/2(7rf) 3 / 2 V 4 ; 



8irt 2 



cos(8t) 



(A39) 
(A40) 



which display a different leading asymptotic behavior compared to the non-equilibrium case after the quench in 
Eqs. (TA"37] and ([MS)) . 
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Appendix B: Global transverse magnetization 



By adopting the same approach and notation as those used for the calculation of the (self-)correlation and (self- 
)rcsponsc function of <5f(i) summarized in App. \K\ we calculate the correlation and linear response functions of the 
total magnetization: 



L 

Its connected correlation function is 

(M*(t 2 )M*(h)) - (M z (t 2 ))(M z (h)) = ^ X>£(iiM*2K(*iK(* 2 ) 



(Bl) 



(B2) 



After a multiplication of this correlation by a factor L — which is required in order for the fluctuations not to vanish 
in the thermodynamic limit — one obtains the symmetric and antisymmetric correlation for L — > oo: 



C™(t 2 ,t 1 )=L 



L 



({M*(t 2 ),M*(h)}) - (M^2))(M Z (M> 



L— voo 



dk 



Re 



Vk{tl)Vk(t2)u* k (h)Uk(t2) 



C M (t2,h) = -{[M'folM'ih)]) L -^8i I — Im ^(tiKfeK(iiK(i 2 ) 



* dk 



(B3) 



As we are primarily concerned with the behavior in the long-time stationary regime, we focus on the stationary 
functions C±(t) = limj^oo C± (t 2 = t + h,ti). The response function R M is related to C* 7 via Eq. ([12]). i.e., 
R (t) = 2i6(t) C_ (t). Inserting the expression for the matrix elements in terms of the Bogoliubov angles, the 
stationary parts of the correlations in Eqs. (|B3I) become: 



n dk 

Gf (t) = 2 / — { (1 - cos 2 A fc ) cos 2 (2# fc ) + [1 - cos 2 (2# fc ) 



cos 2 A fe + 1 



cos(2e fc <) 



R M (t) =i6(t) 



f* dk 

/ _[i -cos 2 (26> fc )]sin(2e^)cosA fc . 
Jo t 



(B4) 
(B5) 



While these expressions are valid for a generic quench we will focus on the critical case, for which e& = t\ =1 = 4 sin(fc/2) 
[Eq. (|34p], cos(2#fc) = sin(fc/2) [Eq. (|A11[) ] and cosAfc can be expressed as in Eq. (|A12[) . In particular, substituting 
into Eqs. (|B4|) and (|B5[) this expression for cos(20fc), renders Eqs. (|64)) and (|65|) which, in turn, can be expressed in 
terms of the integrals in Eqs. (|A16[) . (|A17|) . (|A18I) and of the following ones: 



— sin 2 (fc/2)cos(2e fc t) = - — Ji(8i) 
7r 8 at 

f n dk 

/ — cos 2 A fe cos(2e fe i) = F(8t), 
Jo I" 

dk 

/ — sin 2 (fc/2)cos 2 A fe cos(2e fc t) = -F"{8t), 
Jo 7T 

/" r dk 

/ — sin 2 (fc/2)sin(2e fe f)cosA fe = -E"{8t). 
Jo t 



J (8t) - J 2 (8t) 



(B6) 
(B7) 
(B8) 
(B9) 



We note that the connected symmetric correlation in Eq. (|B4[) contains a constant term, which is attained for t — > oo 
and which can also be expressed in terms of the previous functions [see Eqs. (|B6[) and (|B8[) ] as 



f n dk r i 
C = 2 — (1 -cos 2 A fe )cos 2 (20 fc ) =1 + 2F"(0) 

JO T L J 

(where we used the fact that Jo(0) = 1 and J 2 (0) = 0). Accordingly, one has 

F(8t) + F"(8t), 



C™(t) = C + 
R M (t) = 4 9(t)[E(8t)+E"{8t)] 



(BIO) 



(Bll) 
(B12) 
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with [see Eq. (|B7 



r dfc , , , xx Tsin 2 (/c/2) 

-™<r^ { k,2 )) i + {T _^J [k/2) . (B13) 



The right panel of Fig. [16] shows F(t) as a function of the time r for various values of IV The leading asymptotic 
behavior of the function for t>1 — which is discussed further below in App. IB 11 see Eq. (|B24[) — is indicated by 
dashed lines. A straightforward calculation of F"(0) from Eq. (|B13[) gives [see Eq. (|B10I) ] 



C = (B14) 

(1 + v^) 2 



whereas 



In order to define an effective temperature according to Eq. ()25|) . we will be interested in 



F(0) = — ^— . (B15) 



C¥ (0) - Cf (i = -co) = i + F(0) + F"(0) = ^-=- + — '■= , (B16) 

which follows from Eqs. (|B11|) . (|B10|) . (|B14[) and the fact that Jo, 2 and F vanish for large arguments (as we will 
discuss in App. IB ll below), whereas Jo(0) = 1 and J 2 (0) = 0. Note that in the fully polarized case To = or To = 00, 
both corresponding to T = 1 (see Eqs. (|A13[) and (|A12[1 ), the integrals on the l.h.s. of Eqs. (|B7[) and (|B6[) are identical 
and therefore 

F(t) = [J (r) - J 2 (r)]/2, with F"(r) = -(3/8)J (r) + (l/2)J 2 (r) - (l/8)J 4 (r) - (1 - 3/r 2 )J 2 (r). (B17) 

Taking into account this relation, Eqs. (|B14j) . (|A26|) . (|B11|) . (|B12|) and the properties of the Bessel functions (see, 
e.g., Ref. [64|), one finds 

Cf (t) = J + |jb(8t) + ij 2 (8i) - ij 4 (8t), (B18) 

R M (t) = e(t)[J!(8t) + J 3 (8t)} = e(t)^^-, (B19) 

for a critical quench starting from a fully polarized case To = (or, equivalcntly, Tq — 00). 

In the equilibrium case To = 1(= T), T — > 00 and [see Eq. (|A12|) ] cosAfc(r = l,To = 1) = 1. Accordingly, the 
integral on the l.h.s. of Eq. (|B7|) becomes identical (up to a trivial rescaling of i) to the one on the l.h.s. of Eq. (|A16[) . 
so that 

F(t) = J (r), with F"(t) = -[J (t) - J 2 (r)]/2. (B20) 

Taking into accoun t Eq. (T A29]), o ne fin ds S"(r) = -ffii(r) = i?_ 2 (r) +r- 1 lf_i(r) (see, e.g., 11.4.27 in Ref. @) and 
therefore [see Eqs. (|BTT|) , (|BT2|) . (|BT4|) and 11.4.23 in Ref. [H] 

Cf(t) = Jb(8t) + J2(8i), (B21) 

1 H-!(8ty 



R M {t) = 6{t) 

in the equilibrium critical case at zero temperature (no quench). 



(B22) 



1. Long-time behavior 

In order to determine the long-time behavior of and R M in Eq. (jBllI) and (|B12[) . one has first to determine the 
long-time behavior of the function F introduced in Eqs. (|B7|) and (|B13j) . After the change of variable y = 1 — sin(A;/2), 
F can be cast in the form 

f 1 co S ((l-y)r) 2T (1 - yf 

F(T) = Jo dy Vv -V2- y[ i + (T - m -yry (B23) 
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which has the same form as the real part of the integral analyzed in Eq. (|A35[) . with /(0) = \/2/ir and /'(0) = 
(1 — 8T _1 )/(2v / 27r). Accordingly, the asymptotic expansion of F(t) is given by 

F{T> " W^ sh ' (T + I > + ? W ^ C0S(T + T> + 0( ^ 3) - (B24) 

Taking into account Eqs. (|B24[) , (|A36|i . and the standard expansion of Jo,2( r ) ( scc , e.g., 10.17.3 in Ref. [iij]), one 
eventually finds from Eqs. (jBlljl and (|B12|) : 

Cf (t) = C - ^= sin(J - 8t)^ + G(t- 5 / 2 ), (B25) 

R«(t) = -^=cos(J - 8t)-^ + Or 5/3 )- (B26) 

Note that the case of equilibrium dynamics at zero temperature can be recovered from these expressions by setting 
To = r = 1, i.e., for T _1 = 0. As we discussed after Eq. (|A36|i . the asymptotic behavior of E(t) changes due to the 
additional contribution of a term oc r" 1 which renders the leading decay of R M (t) slower and oc i — 1 , whereas the one 
of C+ 1 is not altered compared to Eq. (|B25[) and — apart from the specific value of C — equal to the non-equilibrium 
one. 



2. Fourier transform of the response and correlation functions 

In order to determine the frequency-dependent effective temperature defined in Eq. (|19[) we calculate here the 
Fourier transforms of and R M in Eqs. (jBllI) and (|B12[) according to the convention (fl~5|) . Consider first C+: 
the functions Jo, J h [see Eq. (|B6|) ] and F [Eqs. (|B7[) and ()B13j) ] in terms of which C+ is expressed involve all an 
integration over the wave- vector k of cos(2efci). As expected from the fact that C¥{t) = C¥{— t), this function is 

even under time reversal. Accordingly, the Fourier transform C+(oj) = dt e lu)i C+ 1 (t) is real and involves integrals 
over momenta of f^^dte lu ' t cos(2ekt) = Tr[S(uj + 2ek) + S(uj — 2efe)], which can be calculated straightforwardly and give 
Eq. (|7Tj) . The response function R M , instead, is given in Eq. ()B12|) in terms of the function E defined in Eqs. (|A18[) 
and (|A22]) . Note that its Fourier transform R M (uj) = f^dt e luJt R M {t) can be also expressed as 

R M (u) = 1 ~ E>(u/8) - where fZ>(w) = / dre^ T S(r) (B27) 



is the Fourier transform of 9(t)E(t) and £"(0) is given by Eq. (IA23[) . Taking into account that E{t) involves an 
integral over momenta of sin(e/ £ r) [scc Eq. (|A22[) ]. with Ek = sin(fc/2), E > (u>) involves an integral over momenta of 



dre l " T sin(£ fe r) = / dr / sin((o; + £fc)T) - sin((o; - £ fc )r) | . cos((o; - e fc )r) - cosQ + g fc )r) 



(B28) 



1111 .Trr/ s r/ M 

-pv — ■ -pv + i-[6(u) - Efc) - d(u) + s k )\, 

2 uj + ek 2 uj — Ek 2 



(see, e.g., Ref. |67| ) where "pv" indicates that the principal value of the subsequent integral over momenta has to be 
considered. Accordingly, from Eq. (|A22|) . 



a M VT f de 

E>{u) = / , ^ —f= = = pv — — - pv- 



7T J y/l^ y/l + (T - l)e 2 \ uj + e uj - e 



+ i- = - VT g(i-H). 
x/T^TJ 2 ^/i + (t - i)w 2 v 1 17 



(B29) 



As expected, due to the fact that E > (uj) is the Fourier transform of the "causal" function 9(t)E(t) which vanishes 
for t < 0, the real and imaginary parts of E > (uj) are connected by a Kramers-Kronig relation [67| . In addition, for 
this specific case, the real part of E > (uj) can be cast in the form 

Rc (uj) = ^ f 1 - d£ 1 (l-oj 2 pv^— ^ , (B30) 
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in which one recognizes the definitions of the complete elliptic integrals of the first and third kind, K e (k) and H(a, k), 
respectively (see 19.2.4, 19.2.7, and 19.2.8 in Rcf. [64]), which we report in Eq. (|A24|) and here for convenience: 



U(a,k) = [ 
Jo 



pv * (B31) 



VT^ 2 VI - ke 2 



(See also the remark after Eq. (|A25[) about our usage of k.) Note that for a > 1 the integrand of II is singular within 
the domain of integration and the principal value has to be considered as a part of the definition. In particular, one 
finds n(a, k) = K e (k) — H(k/a, k) (see 19.6.5 of Ref. [HI), which connects the behavior for a > 1 to the one for a < 1 
(assuming k < 1). Additional properties which will be useful later are: 11(0, k) = K e (k), II(a, 0) = n/(2y/l — a 2 ), 
and K e (Q) = tt/2. According to Eqs. (|A"24l and (lB3"TT) . one can express Eq. (|B30|) as 



Re%M = ^(^ (1 - T ^ II(W " 2 ' 1 - T) fOT | W | > (B32) 
TT |n((l-T)w 2 ,l-T) for M<1, 



in terms of the elliptic integral II(a, k) with a < 1. 

The imaginary part of R M (ui) — which enters into the definition of the effective temperature via Eq. (fT9"]) — can 
be calculated straightforwardly from the imaginary part of E > (u>) above and gives Eq. (1721) . 

In order to define the effective temperature according to Eq. (|25[) one needs to determine R (ui = 0), which is 
related to £>(u; = 0) = ReE>(cj = 0) in Eq. (|B27|) via Eqs. (|B2"9)l and (|B32|) : 

S>(0) - ^^n(0, 1 - T) = ^-KJl - T), (B33) 

7T 7T 

where K e (k) is defined in Eq. (|A24[) (see, e.g., chapter 19 in Ref. [64| and the remarks after Eq. (|A25[) ). Accordingly, 
taking into account also Eqs. (|F327|) and (|A23[) , one eventually finds: 

= 0) = ^^(l-T)- g .(l-T) (B34) 



Appendix C: Short-time expansion of the order parameter correlations 



We aim at calculating C x (t + t , t ) = (o|(t + io)oj(to)) = r (0| a x {t + t )a*{to) |0) r for t <C t <C to an d t given 
in Eq. (|86j) . For the purpose of the present discussion we consider the Schrodinger picture of the dynamics and we 
indicate by |^(to)) = e~ im ° |0) r the state obtained by evolving the ground state of H{^o) with H = H(r) up to the 
time to- Expanding the evolution operator for small £ one obtains: 



C*(t + to,to) = r o {0\e iA{t+to) ct?e- i6t a*e- i6t ° |0> r „ = (V»(*o)| e^e"****? |^(t )) 

1 



(V(to)| i + *t[ff,5?]s? + -* i 



oj + 0(^) |V(*o)>, 



(CI) 



where 1 is the identity matrix. The commutators in the previous expression can be calculated by taking into account 
that Pauli matrices at different sites commute, whereas ot.a^ = S a i,l + ie 



tensor with e xy 



&j, with e 

1. Accordingly, setting J = 1 in the definition of H in Eq. 



-abc 



abc the completely antisymmetric 
26]) . the following commutation 



relations hold: [H,a?]df = -2Ta^ and {H,a*},H a* = -4T 2 + 4T a^a^a^ + af +1 ). We note that the last term 

of the latter expression is an anti-hcrmitcan operator and therefore its expectation value is imaginary, contributing 
only with a subleading, 0(t 2 ) term to the response function cx ImC x , which receives a contribution 0(t) from the 
first commutator calculated above. In addition, one can actually verify that such a term vanishes in the stationary 



regime. In fact, by using Eq. (|27[) one finds oj <7j +1 
and ([201) 



c]ct +1 + ctcj +1 +h.c. = (cj — Cj)(cj , x +Cj+i) and from Eqs. ([H 



k,k' 



3 -t(fc+fc')i+i*atftt _ p -i(fc-fe'b'+»*»t 



-ft 1 * 



C fc C fc' 



(C2) 



43 



and analogously 

~ OjOiof+i = ( £ j + cj)(c]+i + 

1 r i (C3) 

fe,fe' 

Upon adding Eqs. (jC2[) and (|C3|) and rearranging the terms by using the canonical anticommutation relations of the 
fermions one obtains: 

= jY^{( iaink ) e ~ ilk+k ' )j ^k' +e _-i(fc+fe ' w (cosfe'-cosfc)44' +e i(fe+fe ' w (-isinfc)c fe c fe /|. 
fe.fe' 

Its expectation value is therefore given by 

- (V(to)| +of+i) |V(*o)> = ^^(sin^Re^^oK^o)] = * ]T sin fc sin A fc sin(2e fc i ), (C5) 

fe>0 fe>0 

where we used the fact that (V(*o)l44' l"0(*o)) = r o (0| 4(*o)4'(*o) l°)r = -v k (to)u* k (t )5 k >- k , (VK*o)l4 e fc' lV>(*o)) = 
|wfc(to)| 2 4',fe, and (V>(*o)l c fc Cfe' |-0(*o)> = Wfc(*o)Vfc(*o)*fc,-fc', which follow from Eqs. ([37]), ([55]) and (fM} (in which we 
omit the superscripts of u k , v k ), with e k = e£ and A& = 2(0^ — 6>£ ). The expression on the r.h.s. of Eq. (|C5[) can be 
shown to vanish in the stationary regime to oo. Therefore, the expressions above give for t, to > 0: 

C*(i + i ,*o) = ReC x (t + t ,t ) = 1- ^AT 2 +0(t 3 ), 

(C6) 

R*(t + t ,t ) = -2lmC x (t + t ,t ) = 4Tt (V(to)|<x? |^(*o)> + 0(« 2 ), 

which, in the stationary limit to — > oo, yield Eqs. (|82[) and (|83[) where, in indicating the corrections, we used the fact 
that one expects the correlation function in the stationary regime to be an even function of time. 



Appendix D: Energy-based effective temperature Tjj(r, To) in the limit of shallow critical quenches 

The effective temperature T® s based on the energy of the system is obtained according to Eq. (|45)l , i.e., by equating 
the thermal average of the energy at a certain temperature Tj^, see Eq. (|44|). with the expectation value in Eq. (|40|) 
of the energy after the quench from To to Y. Here we are interested in a quench of the Ising model in a transverse 
field r with critical final value r = 1. As discussed in Sec. IIII CI the expectation value of the energy after the quench 
is given by Eq. (|44|) : 

fx Ah 

r o (0|#(r)|0) ro = - / — e k (T = l)cosA fc (r = l,r ) 

r- ( D1 ) 

4^T EJl - T) - K e (l - T) 
= -2E' = -— ^ '-, 

7T T — 1 

where we used the fact that, at the critical point, efc(T = 1) = 4| sin(fc/2)| and wc took into account the definitions of 
the functions cosAfe(r = 1,Tq) and E'(0) given in Eqs. (|A12|) . (|A22|) . and (|A23|) in order to express the result as a 
function of T, see Eq. (|54|) . For Fo — > 1, T — > oo and therefore one can expand the expression above: 

ro (0|^(r = l)|0) ro = -1 + ' 3 + 41n T 2 + lnT + G(T- 2 ,T- 2 lnT). (D2) 

This expansion has to be compared with the thermal average of H(T = 1) at a temperature /3 _1 = Tj^, which is 
expected to vanish for To — > 1. The thermal average given by Eq. (|44p can be conveniently cast in the form 

dk 

(HQ? = l)) T= p-i = - / — e k (T = 1) tanh ((5e k (T = l)/2) 



d dk 

-2 — / — lncosh(2p , sin(fc/2)) (D3) 
dp Jo 2?r 

7r dp J 2n \ ) 
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In the limit (3 — > oo we are interested in, the last integral is dominated by small values of k, and therefore one can 
approximate sin(/s/2) ~ k/2 in the integrand: 

where the extension of the integral from [0,7r] to [0,oo] introduces only exponentially small corrections. Accordingly, 

(^(r-DW/- 00 -^^, (D5) 

and, by comparing Eqs. (|D2[) and (|D5|) , one finds: 

oE _^ \ \ ^ ^ (D6) 

° ffl r e f(r=l,r ^-l) ~ 2 A /6(lnT + 41n2-3)" 1 ' 

This asymptotic expression is used in Fig. [I5]in order to plot the behavior of Tj^r = 1, To) for To > 0.9. In addition, 
this result can be compared with the inverse temperature obtained by studying the zero-frequency limit of the FDT 
of the global transverse magnetization M (see Sec. II V Bl and in particular Eq. ((74")) ) 



which shows that T c ^/T c f -> for T a -> T = 1. 



£ = ^ = ^ (D7) 



[1] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, Rev. Mod. Phys. 83, 863 (2011); J. Dziarmaga, Adv. in Phys. 
59, 1063 (2010); A. Dutta, U. Divakaran, D. Sen, B. K. Chakrabarti, T. F. Rosenbaum, and G. Aeppli. larXiv : 1012 . 0653 
(2010). 

[2] M. Greiner, O. Mandel, T. W. Hansen, and I. Bloch, Nature 419, 51 (2002). 

[3] S. Trotzky, Y.-A. Chen, A. Flesch, I. P. McCulloch, U. Schollwock, J. Eisert, and I. Bloch, Nature Phys. 8, 325 (2012). 
[4] M. Cheneau, P. Barmettler, D. Poletti, M. Endres, P. Schauss, T. Fukuhara, C. Gross, I. Bloch, C. Kollath, and S. Kuhr, 

Nature 481, 484 (2012). 
[5] L. Foini, L. F. Cugliandolo, and A. Gambassi, Phys. Rev. B 84, 212404 (2011). 
[6] G. Gallavotti, Statistical mechanics: A short treatise, (Springer Verlag, Berlin, 1999). 
[7] J. von Neumann, Z. Phys. A 57, 30 (1929). 

[8] S. Goldstein, J. L. Lebowitz, R. Tumulka, and N. Zanghi, Eur. Phys. J. H 35, 173 (2010); S. Goldstein, J. L. Lebowitz, C. 
Mastrodonato, R. Tumulka, and N. Zanghi, Phys. Rev. E 81, 011109 (2010). 

[9] H. Tasaki, Phys. Rev. Lett. 80, 1373 (1998) : laxXiv : 1003 ■ 54231 (2010) . 
[10] M. Srednicki, Phys. Rev. E 50, 888 (1994). 
[11] J. M. Deutsch, Phys. Rev. A 43, 2046 (1991). 

[12] L. F. Cugliandolo, J. Kurchan, and L. Peliti, Phys. Rev. E 55, 3898 (1997); L. F. Cugliandolo and J. Kurchan, J. Phys. 

Soc. Japan 69, 247 (2000). 
[13] A. Crisanti and F. Ritort, J. Phys. A 36, R181 (2003); L. Leuzzi, J. Non-Cryst. Sol. 355, 686 (2009). 
[14] L. F. Cugliandolo, J. Phys. A 44, 483001 (2011). 

[15] A notable exception to this scenario is provided by classical systems during critical coarsening, for which the effective 
temperature that can be extracted from the dynamics via FDRs turns out to depend on the choice of the observable 
beyond the mean-field approximation. See P. Calabrese and A. Gambassi, J. Stat. Mech. P07013 (2004); J. Stat. Mech. 
P01001 (2007). 

[16] D. Karevski, [arXiv: cond-mat/06113"27l (2006). 

[17] M. Rigol, V. Dunjko, V. Yurovsky, and M. Olshanii, Phys. Rev. Lett. 98, 050405 (2007); M. Rigol, A. Muramatsu, and 

M. Olshanii, Phys. Rev. A 74, 053616 (2006). 
[18] M. A. Cazalilla, Phys. Rev. Lett. 97, 156403 (2006); A. Iucci and M. A. Cazalilla, Phys. Rev. A 80, 063619 (2009). 
[19] M. A. Cazalilla, A. Iucci, and M.-C. Chung, Phys. Rev. E 85, 011133 (2012). 
[20] D. Fioretto and G. Mussardo, New J. Phys. 12, 055015 (2010). 
[21] P. Calabrese, F. Essler, and M. Fagotti, Phys. Rev. Lett. 106, 227203 (2011). 
[22] B. Blass, H. Rieger, and F. Igloi. larXiv : 1205 ■ 33031 (2012) . 

[23] P. Calabrese, F. H. L. Essler, and M. Fagotti, J. Stat. Mech. (2012) P07016; J. Stat. Mech. (2012) P07022. 

[24] M. Rigol and M. Srednicki, Phys. Rev. Lett. 108, 110601 (2012). 

[25] C. Kollath, A. Laeuchli, and E. Altman, Phys. Rev. Lett. 98, 180601 (2007). 



45 



G. Biroli, C. Kollath, and A. Laeuchli, Phys. Rev. Lett. 105, 250401 (2010). 

T. Caneva, E. Canovi, D. Rossini, G. E. Santoro, and A. Silva, J. Stat. Mech. P07015 (2011). 

M. C. Banuls, J. I. Cirac, and M. B. Hastings, Phys. Rev. Lett. 106, 050405 (2011). 

G. Carleo, F. Becca, M. Schiro, and M. Fabrizio, Sci. Rep. 2, 243 (2012). 

D. Rossini, A. Silva, G. Mussardo, and G. E. Santoro, Phys. Rev. Lett. 102, 127204 (2009). 

D. Rossini, S. Suzuki, G. Mussardo, G. E. Santoro, and A. Silva, Phys. Rev. B 82, 144302 (2010). 
R. Barnett, A. Polkovnikov, and M. Vengalattore, Phys. Rev. A 84, 023606 (2011). 

M. Eckstein, M. Kollar, and P. Werner, Phys. Rev. Lett. 103, 056403 (2009). 
S. Sotiriadis and J. Cardy, Phys. Rev. B 81, 134305 (2010). 

M. Schiro and M. Fabrizio, Phys. Rev. Lett. 105, 076401 (2010); Phys. Rev. B 83, 165105 (2011). 
P. Calabrese and J. Cardy, Phys. Rev. Lett. 96, 136801 (2006); J. Stat. Mech. P06008 (2007). 
S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, 1999). 
J. Cardy, to appear (2012). 

Super-cooled liquids are systems in equilibrium (metastable with respect to the crystal) with the thermal bath in which 
the decay of the correlation functions is rather complex, certainly non-exponential and for which the identification of a 
single coherence time is delicate. 

C. Godreche and J.-M. Luck, J. Phys. A 33, 1151 (2000). 

P. Calabrese and A. Gambassi, J. Phys. A 38, R133 (2005); F. Corberi, E. Lippiello, and M. Zannetti, J. Stat. Mech. 
P07002 (2007). 

M. Henkel and M. Pleimling, Nonequilibrium Phase Transitions Vol. 2 - Ageing and Dynamical Scaling far from Equilibrium 
(Springer, Heidelberg, 2010) 

L. F. Cugliandolo and G. S. Lozano, Phys. Rev. Lett. 80, 4979 (1998); Phys. Rev. B 59, 915 (1999). 

F. Igloi and H. Rieger, Phys. Rev. Lett. 85, 3233 (2000). 

F. Igloi and H. Rieger, Phys. Rev. Lett. 106, 035701 (2011). 

F. Igloi and H. Rieger, Phys. Rev. B 84, 165117 (2011). 

J. Marino and A. Silva, [arXiv : 1203 ■ 21081 (2012') 

A. Silva, Phys. Rev. Lett. 101, 120603 (2008). 

A. Gambassi and A. Silva, larXTv: 1106 ■ 267T1 (2011). 

V. Mukherjee, U. Divakaran, A. Dutta, and D. Sen, Phys. Rev. B 76, 174303 (2007). 

A. Caso, L. Arrachea, and G. S. Lozano, Phys. Rev. B 83, 165419 (2011); Phys. Rev. B 81, 041301(R) (2010); 
larXiv: 1111 .19881 (2011). 

S. Sotiriadis, P. Calabrese, and J. Cardy, EPL 87, 20002 (2009). 

R. Kubo, M. Toda, and N. Hashitume, Nonequilibrium Statistical Mechanics, 2nd ed. (Springer Verlag, Berlin, 1991). 

B. Sciolla and G. Biroli, Phys. Rev. Lett. 105, 220401 (2010); J. Stat. Mech. P11003 (2011). 

E. Altman and A. Auerbach, Phys. Rev. Lett. 89, 250404 (2002). 

R. A. Barankov and L. S. Levitov, Phys. Rev. Lett. 96, 230403 (2006); E. A. Yuzbashyan and M. Dzero, Phys. Rev. Lett. 
96, 230404 (2006); E. A. Yuzbashyan, O. Tsyplyatyev, and B. L. Altshuler, Phys. Rev. Lett. 96, 097005 (2006). 
A. Gambassi and P. Calabrese, EPL 95, 66007 (2011). 

A. Mitra and T. Giamarchi, Phys. Rev. Lett. 107, 150602 (2011); Phys. Rev. B 85, 075117 (2012); J. Lancaster, T. 
Giamarchi, and A. Mitra, Phys. Rev. B 84, 075143 (2011). 

I. Lesanovsky, B. Olmos, and J. P. Garrahan, Phys. Rev. Lett. 105, 100603 (2010). 

B. M. McCoy, E. Barouch, and D. B. Abraham, Phys. Rev. A 4, 2331 (1971). 
E. Lieb, T. Schultz, and D. Mattis, Ann. Phys. 16, 407 (1961). 

E. T. Jaynes, Phys. Rev. 106, 620 (1957); Phys. Rev. 108, 171 (1957). 
E. Barouch, B. M. McCoy, and M. Dresden, Phys. Rev. A 2, 1075 (1970). 

NIST Handbook of Mathematical Functions, edited by F.W.J. Olver, D.W. Lozier, R.F . Boisvert, and C.W. Clark (Cam- 
bridge University Press, 2010). See also the on-line version http://dlmf.nist.gov/ which we refer to when quoting 
formulas and definitions. 
T. Niemeijer, Physica 36, 377 (1967). 
P. Mazur, Physica 43, 533 (1969). 

W. Appel, Mathematics for Physics and Physicists (Princeton University Press, 2007). 

D. Schuricht and F. Essler, J. Stat. Mech. P04017 (2012). 

J. H. H. Perk and H. Au-Yang, J. Stat. Phys. 135, 599 (2009). 

L. F. Cugliandolo, D. R. Grempel, G. S. Lozano, H. Lozza, and C. A. da Silva Santos, Phys. Rev. B 66, 014444 (2002). 
For some particular values of V and Fq, we verified numerically that the the long-time exponential decay rate of C+(t) 
and R x (t) in the stationary regime is well-fitted by Eq. (|86[) also for non-critical quenches r / 1. This suggest that the 
conclusion we draw here for r = 1 might extend also to the case r / 1. 

P. Deift and X. Zhou, in Singular limits of dispersive waves (Lyon, 1991), p. 183, NATO Adv. Sci. Inst. Ser. B Phys. 320 
(Plenum, New York, 1994). 

S. Sachdev and A. P. Young, Phys. Rev. Lett. 78, 2220 (1997). 

D. Rossini, S. Suzuki, and G. E. Santoro. larXiv : 0910 . $0551 (2009) . 



